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The general problem chosen for consideration concerns the 
establishment of a payload in orbit about the planet Mars. Empha- 
Sis in this study has been placed upon a comparison of the various 
paths available for a free-fall, space stabilized rocket, the re- 
quirements necessary to achieve such paths, and the selection 
of the optimum path. The optirnum path has been compared with 
other trajectories and conclusions have been drawn as to the 
requirements for deviations from the optimum path. The study 
is primarily concerned with the midcourse portion of the trip, 
which has been defined by the authors as that part of the overall 
trajectory extending from a desired circular orbit about Earth 
to a desired circular orbit about Mars. A firm basis for the 
se:ection of the orbit has been established. Treatment has also 
been given in this report to the cortinuous determination of 


present position 1n space and ‘he initiation of corrective action 


Cao 


when deviation from a desired path exists. The study also 
includes generalized instrumentation and guidance requirements. 
Extensive quantitative inforrnation has heen compiled which will 
serve as a basis for the determination of instrumentation re- 
quirements as well as creata a proper background for an 


appreciation of flight within the solar system. 
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OBJECT 


The object of this paper is to investigate the mid-course 


guidance requirements to place a payload in orbit about Mars. 


Char i) iva 


INTRODUCTION 


The authors have undertaken this study in order to bring 
together, in proper context, some considerations which may be 
anticipated in interplanetary flight. A journey to the planet Mars 
has been chosen as the basis for the presentation. However, 
many of the ideas and methods set forth here can be considered 
applicable to flight to other bodies within the solar system. It 
is felt that no loss of generality will result if attention is focused 


upon this particular choice. 


This paper has evolved as a study of several problems as- 
sociated with interplanetary flight. The objective to be pursued 
is concerned with placing a payload in orbit about the planet Mars. 
Departure is assumed to originate from orbit about Earth. The 
attainment or possibility of flight is not questioned. Rather, 
this is a broad study of some preferred methods, minimum re- 
quirements, and ideas concerned with fulfillment of the general 


Oblective. 


Recent developments have justified studies of interplanetary 
flight. It is realized, however. that the topic is hopefully opt: - 
mistic in many respects. For this reason the study is neces- 
sarily broad and is not designed to give an intricate, detailed 


account of methodology 


A flight to the moon will logically precede planetary flight, 
but its close proximity causes this problem to become quite 
specialized as contrasted to the Mars flight. Piloting a boat in 
an inland waterway as compared to navigating a ship on the open 


sea is an appropriate analogy. The two problems are independent 


in a large sense. : 
This study has been preceded by an extensive search through 

the literature for information regarding this comparatively new 

field. Treatments have generally varied from the popularized 

texts to the highly specialized mathematical devedopments. 

There appears to be a paucity of practical approaches to the 

problem which can be interpreted by the systems deSigner or 


engineer. 


The authors generally envisage an unmanned vehicle with 
small payload, unpowered except for terminal and corrective 
thrust requirements. It is believed instrumentation can be 
handled with components whose performances do not differ too 


radically from slight extensions of present day development. 


A practical simplified view has been taken as a rule due 
to the lack of an available high speed digital computer. Such 
a device would of necessity be required if accurate trajectory 
Studies were inorder. However, it is felt that little departure 
from reality exists under the assumptions imposed. In fact, 


the simplifications add to the comprehension of the material. 


Solutions are sought in answer to the following questions: 
(1) What are the major influences of the solar system 


upon the rocket in flight? 
(2) What simplifying assumptions can be made? 
(3) What routes may be taken? 


(4) What criteria are available for the selection of 


one trajectory or route in preference to another ? 


(5) How may position in space be continuously 


determined for navigational purposes? 


(6) How may the rocket be assured of arriving in the 


target area if errors do occur? 


Ze 


I 


CHAPTER 2 


THE SOLAR SYSTEM 


A summary of currently accepted data concerning the 
universe, which is pertinant to this study, is needed for formu- 
lation of assumptions and as basis for computation. This is 
presented below as oetesion and in tabular form and figures, 


and has been derived from references” 49,44, 33, 94, 50" 


The Sun and its Planets 


The solar system consists of the sun, nine conspicuous 
planets, and numerous smaller planets. The conspicuous planets 
are placed in two groups. The major planets, Jupiter, Saturn, 
Uranus, and Neptune, are large, have low density, and are 
surrounded by high atmospheres. The terrestrial or minor 
planets, Mercury, Venus, Earth, Mars, and Pluto are similar 
to the earth in many respects. The so called navigational planets 
are Venus, Mars, Jupiter, and Saturn, deriving this classi- 


fication from their nautical use for position determination. ” 


For the period of 1956-1970, the angular motion of the 
major planets in the solar system will be as shown in Fig. 2.1 
The minor planets, Mars, Earth, and Venus, are shown in their 
positions as of 1 January 1970, and in their respective orbits, 


drawn as circles. 


A generalization which serves to illustrate the area occupied 
by the planets and sun is to picture a piece of 1/4 inch plywood cut 
in a circle of 5 feet diameter. This would contain all the planets 
in their orbits around the sun (except Pluto) if drawn to a scale 
such that Neptune's major planetary axis was 9 feet. 


* Superscript numbers refer to bibliography. 
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Fig.2.1 Angular motion of major planets from 1965 to 1970 
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solar Radiation and Temperature in Space 


The solar constant of sun's radiation is 1. 938 calories per 


minute per Square centimeter, +.05 — This value is equiva- 


lent to 1.35 x 10° ergs per second per centimeter square or 
otal x ioe ergs per second per inch Square at a distance of 1 AU 
from the sun. This incident energy varies as the reciprocal of 


the range squared, as shown in Fig. 2.2 


The maximum temperature of the surface of a black body 
planet with no atmosphere varies as the reciprocal of the square 
root of the distance of the body from the sun, as shown in fig.2. 3. 


At a distance of 1 AU from the sun this temperature is 392°K. 


Approximate Orbital Paths of Planets 


The path of a planet or body in orbit about the sun can be 
completely described by seven elements, in so far as the body 


depends only upon the sun's attractive force. These are: 
1. Semi major axis (a) 
2. Eccentricity (e) 
oleae latte 1 vaya (Gi) 
4. Longitude of ascending node (23) 
5. Longitude of perihelion (7) 
6. Period (T) 
7. Epoch 


These are shown in Fig. 2.4 and 2.5. 


Kepler's Three Laws of Planetary Motion are also appli- 


cable and are: 


1. The planets move in ellipses, with sun at one of 
the focuses 
2. The radius vector sweeps out equal areas in equal 


time 
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Fig. 2.2 Incident solar energy as a function of distance from Sun 
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Fig. 2.4 Basic ellipse 
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Fig. 2.9 Basic ellipse oriented to reference plane 


3. The square of the time of revolution is proportional 


to the cube of the mean distance. 


2.4 Perturbations of Planetary Orbits 


The actual motion of the planets may best be described as 


follows: 


1. The orbits of the planets are not fixed but gradually 
change in eccentricity, inclination, (etc.). These 
changes are classified as secular perturbations; 
however, these perturbations are so slow that during 
any one revolution the orbit is altered very little and 
during the time of any given single interplanetary 


flight, can be neglected. 


2. A planet in its motion does not exactly follow 
Kepler's Laws but is seen to vary in three dimensions 
relative to its ideal Keplerian Ellipse. Such changes 
are very small and are called periodic perturbation. 
These values expressed in heliocentric longitude A) 


are. 


Mercury 15" 


Venus 30'! 
yas tla ih 
Mars ye 


Jupiter 30! 
Saturn 10! 
Uranus 60! 


Neptune 30! 


In general the larger the terms the longer the period of 


periodic perturbation. For example Jupiter and Saturn exhibit 


periods of periodic perturbations of about 913 years. 


oO 
Co 


The Invariable Plane 


The most stable property of the solar system is the so 
called invariable plane. This is defined by the condition that the 
total angular momentum of the solar system is zero about any 
axis inthis plane. The invariable plane can be considered as 
the fundamental reference plane of the solar system and is 
inclined 1° 38' 7" to the ecliptic reference plane of 1900 with 
ascending node in longitude 106° 42' measured fromY . This 
inclination is decreasing at the rate of .144" per year and will 
reach 47' in about 20, 000 years at which time it will begin to 


increase again. 


Cislunar Motion 


Cislunar (earth-moon) motion in the plane of the ecliptic is 
that of two unequal revolving masses in a central force field, 
whose center of mass follows a slightly distorted Keplerian 
ellipse of eccentricity (e)=.01674. This cislunar center of mass is 
approximately 2880 miles from the earth center in the direction 


of the moon. SeeFig. 2.6 and 2.7. 


The moon's path around the earth has an eccentricity .055 
and a mean inclination to the ecliptic of 5° 9'. The average 
distance from earth to moon is 238, 857 miles (60.207 earth radii) 
with maximum distance of 252, 710 miles and minimum distance of 
221,463 miles. : 


If we represent the cislunar orbit by a circle of radius 100 
feet, the moon will deviate from this circle approximately 3 inches 
on either side in one transit around the earth. During this time 
the geocenter will deviate approximately .04 inches either side 
of the orbit. 


Solar Motion 


Campbell and Moore in 1920 found the solar system to be 


moving at about 12.2 miles per second wiih respect to the nearer 
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The Invariable Plane 
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that of two unequal revolving masses in a central force field, 
whose center of mass follows a slightly distorted Keplerian 
ellipse of eccentricity (e)=.01674. This cislunar center of mass is 
approximately 2880 miles from the earth center in the direction 


of the moon. SeeFig. 2.6 and 2.7. 


The moon's path around the earth has an eccentricity .055 
and a mean inclination to the ecliptic of 5° 9'. The average 
distance from earth to moon is 238, 857 miles (60.207 earth radii) 
with maximum distance of 252, 710 miles and minimum distance of 
221,463 miles. 
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If we represent the cislunar orbit by a circle of radius 100 
feet, the moon will deviate from this circle approximately 3 inches 
on either side in one transit around the earth. During this time 
the geocenter will deviate approximately .04 inches either side 
of the orbit. 
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few hundred stars toward the constellation Hercules. Jay, 
sanford, Stromberg, and Hubble, among others, estimate that 
the solar system and the neighboring stars are moving at about 


170 miles per second toward the constellation Cygnus. 


Conclusions reached by Oort place the solar system about 
32,000 light years from the center of the galaxy which is in 
Sagittarius and that this position is about two thirds of the dis- 
tance from center to circumference. At 170 miles per second 
the periodof galactic revolution of our sun is approximately 
2,200,000 years, indicating that during the time covered by any 
given flight within the solar system, the sun may be considered 


as a rather stable reference point. 


Fixed Stars 


The nearest star, Alpha Centauri, is at a distance from our 
sun of 4.3 light years or 2.55 x 1G. miles. It has also been de- 
termined that no star exhibits a parallax of more than .76 seconds 
and only 200 stars have proper motion greater than one tenth of 
"Barnard's Runaway Star" which traverses 10.25 seconds of arc 
per year. For the purposes of this paper a ''Fixed" star, as the 
term is used for establishment of celestial inertial reference, 
will be considered to be a star exhibiting a parallax and a proper 
motion of less than one quarter second of arc.Ofthe stars checked 
to date only about .01% do not fall within this Classification of 


"Rixed'' stars. 


Gravitational Effects 


Newton's general Law of Gravitation states that, ‘Any two 
particles of matter attract each other with a force proportional 
to the product of their masses and inversly proportional to the 


N94 or f= KL 
R 


the so called Gaussian gravitational constant equal to 





square of the distance between them where K is 


3 


6.67 x oe ——— From this the acceleration of gravity 
gM(m) sec 
is then = — : 
R 


Values of the acceleration of gravity for the various solar 
bodies are tabulated in Table 2.2 with the mean surface 
gravity for earth equaltol. Fig. 2.8 presents ''g' versus the 
radial distance from the earth, mon and Mars and the variation 
in the effect of the sun's gravitational attraction over the distance 


between the orbits of Earth and Mars. 


Herrick has pointed oe that Gauss was able to determine 
mathematically the value of ke for heliocentric orbits with dis- 
tance in AU, accurate to nine significant figures but neither the 
laboratory value of gravitational acceleration nor laboratory units 
of measurements at present will permit any greater accuracy than 


3 or 4 significant figures. 


An interesting discussion of Gravitational and tidal effects 
which has direct bearing on any travel between Earth and Mars 
has been presented by Schaub”? . Schaub states that the inner 
most moon of Mars, Phobos, is within Rocke!s eee This 
means that since Phobos is assumed to have the same density as 
Mars and is traveling in an orbit of less than 2.44 times Mars 
radius (see Table 2. 2 ) its surface is now free of all loose 
material. Eventually, due to the bombardment by rays from 
space, more and more of this moon will be torn away forming a 
ring of dust particles around Mars. ‘This implies that any 
terminal orbit about Mars should stay well clear of Phobos less it 
hasten this disintegration possibly destroying the vehicle in the 


process. 
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RADIAL DISTANCE FROM } MARS IN MILES 
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Fig.2.8 Gravitational field effects 
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CHAT TE ihe) 


THE PROBLEM SIMPLIFIED 


Introduction 


The Law of Universal Gravitation as enunciated by Newton 
states ‘Every particle of matter in the universe attracts every 
other particle with a force which acts in a line joining them, and 
whose intensity varies as the product of the masses and inversely 


as the squares of their distances apart. "' 


4 


Moulton** presents an excellent discussion of this law and 
states, ‘It will be observed that the law of gravitation involves 
considerably more than can be derived from Kepler's Laws of 
Motion; and it was by a master stroke of genius that Newton 
grasped it in its immense generality and stated it so exactly that 
it has not been necessary to change a single syllable in over two 
hundred years. When contemplated in its entirety it is one of the 


prandest conceptions in the physical sciences." 


It is evident that the various bodies comprising the Solar 
System will exert certain forces upon a rocket in flight among 
them. Because of the extremely large masses and distances 
involved we will omit consideration of the minor bodies, such as 
asteroids and meteor swarms, and show later by inference that 
this omission is justifiable. Considering then only the Sun, Moon, 
rocket, and nine planets, we are faced with a twelve-body problem. 
It would appear expedient to simplify this problem if possible 
Since twelve bodies require thirty-six second order differential 
equations(thirty-three, if referred to the Sun) to adequately de- 


scribe their motions. 


me 


Differential Equation of Motion of the Rocket 


In order to reduce the complexity of the problem let us first 
examine the terms of the equation of motion of the rocket and 


consider their magnitudes. With reference to.Fig. 3.1 define 


I as an inertial reference frame in the convention of Wrigley” 


and 0 as a frame within the Solar System non-rotating with respect 
to I and with axes parallel to those of I. Note that at this point 0 
is restricted only to lie somewhere in the vicinity of the Solar 


System. A represents the position of the rocket, a variable. 


We may write Newton's Second Law in vector notation as 


follows; 


Bee 


[Ryaly = yey fe f, (3-1) 


where [ Rialy is the second derivative of the vector Ria with 


respect to the frame I and the other terms are previously defined. 


nen; 
ora + Ger + Ge + Go, + Gon + Gorn 

+ Gan + Gpa + Ga, + Sys tan (3-2) 
where G, A is defined in Derivation Summary}. 


If there is a particle at point 0 of mass mM then; 


mo(R 1, = ™, (Geo Ge te i a [So ee es: ) 
and the masses may be cancelled. 


By geometry 


[Ryaly = ([Roaly + (Ryo) y (3-4) 
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(a) INFERIOR BODY (b) SUPERIOR BODY 


Fig. 3.2 Planetary configurations for Table 3.1 
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Hence by substitution 


(Rig) ‘ [Ra r= fp + fing ae 
or 
Cont Ca C7 eee ls ley. 
Gey Gaye er, + . +: (3-6) 


Hence in full 
[Ro al t= fae + Gon ~ Gao) + Gon - Gao) + Gea - Seo) 
HGesn ~ Geo) + Gon ~ Go.) Gg - F go) 
HGap, - Gao) + (Gpa - Gao) + (GS, - Go) 
HGya - Gyo) + (Gp, - Feo) (3-7) 


We shall now consider applying this derivation (frame 0) to 


various members of the Solar System. 


Inertial Points Within the Solar System 


Looking first at the motions of the Sun, or the Solar System 
as a unit, it appears from a statistical analysis of the motions 
of the nearer" stars that the Solar System possesses a velocity 
of approximately 12.2 miles per second relative to these "nearer" 
stars. The nearer" stars range from 4 to 300 light-years 
distant from the Sun. Further, the Sun is thought to occupy a 
position approximately 20, 000 light-years from the center of our 
Galaxy (The Milky Way) with a period of revolution about the 
galactic center on the order of 2,200,000 years. (See Section 2.7 
for further details on the Solar System). In addition our Galaxy 
is a member of the so-called Local Group consisting of 17 galax- 


ies, clusters, and clouds including the Andromeda Galaxy, The 


Large and Small Magellanic Clouds, and several small clusters. 

It is thought that these 17 members rotate about a common center 
of mass located somewhere between our galaxy and the Andromeda 
Galaxy. Lastly, current theories on the ''Red Shift" of the optical 
spectra of distant galaxies conclude that these galaxies are re- 
ceding from the Solar System at such tremendous relative ve- 
locities that certain effects attributed to the Theory of Relativity 
must surely be present. It is not the purpose of this paper to 
argue any of the current theories or to attempt any analysis in- 
cluding Relativity Effects. These motions to which the Solar 
system is subject are mentioned solely as a basis for the following 


statements: 


The Sun is probably undergoing motions involving intricate and 
complex accelerations; however, the vast macrocosmic scale of 
distances involved is such that during the time of flight of our 
rocket the effect of these ill-defined accelerations is neither ap- 
preciable nor measurable. This then permits the statement, 
"Newtonian space is that unaccelerated with respect to the 'fixed' 
stars’, i particularly since it is not known that any stars are 
fixed, but only that their motions are very small during recorded 


history. 


The foregoing definition of Newtonian space is preferred by 
the authors since it appears at least semantically more defi- 
nitive than the statement ''Newtonian space is that space wherein 
Newton's Laws are valid.""+ We say then "The possible motions 
of the center of the Solar System, however complex, within our 
detective ability are not known to interfere measurably with the 
information obtained while ignoring these motions and effects 


thereof." 


Simplified Differential Equations of Rocket Motion 


On the foregoing basis we allow ourselves to define the center 
of (1) the Sun, (2) the Earth, and (3) Mars each in turn as the lo- 


ey) 
Ut 


cation of our inertial reference frame 0. Equation (3-7) may 


then be written successively; 


[Roaly = tf ai Goat (Gon -GadttGg an Geo ) ; o 8 ioe 


[Rol = a + Gon + (Gon = Goo) + (Cran < Gad) pa 
3-10 


Table 3.1 gives values in ft/sec’ of the maximum values 
the above terms in parenthesis will hold, using mean values for 
the orbital radii of the bodies. Use of these mean radii is con- 
sidered justifiable since only rarely will the bodies be aligned 
such that these maximum values actually occur. Figure 3.2 


clarifies the manner in which these values do occur. 


By inspection of Table 3.1; if we establish that all acceler- 
ations less than mw” ft/sec” are to be ignored, we may eliminate 
from future consideration most of the parenthetic terms in the 


foregoing equations. We then have; 
UR oalr = fn * Gen HGea ~ Geo) + Fea > Seo) 
+ (Gon- Go ) (3-11) 


[Realy = ff + Gan HGoa - Gog) HIG ¢ A - Gee ) 


HG ot, - Gof) (3-12) 
(Roly = fae + Soy t Gog Gad) + (Gg n- Ge) 
HGar- Geog) (3-13) 


We have now reduced the problem to one of five bodies rather 


than twelve. We will investigate the possibility of reducing it even 
further. 


ELS) 
oy) 


ee 


sae 
en 
1.1] 7 240,000 miles 
vf 1" 
9-2| = 363, 000 
a " 
3-3| = 953, 000 
= 7 iT 
4.4| = 753, 000 





Fig.3.3 Illustration of Rocket-Moon Separation 


The Lunar Gravity Field Effect 


Referring to the Hohmann Ellipse“ it may be shown that in 
the first 15.2 days of rocket flight along this trajectory the radius 
vector of the rocket measured from the center of the Earth 
increases by approximately 600,000 miles. During this period 
the rocket is in almost linear flight, hence the magnitude of the 
Lunar gravity field is easily estimated and is shown in Fig. 3.3. 
At the end of the 15.2 days the value of the gravity field has de- 


creased to less than 1x 107° ft/sec* and is continually decreasing. 


As a result of these rudimentary estimations it is evident 
that the effect of the lunar gravity 11eld need be Considenecsen, 


during the initial portion of the trajectory. The example selected 


above is unique in that it has been shown to require the least possi- 
kle rocket velocity hence the rocket would be subjected to Lunar 
influence for a longer time than in any other trajectory. The fore- 
going estimates then are an upper bound to the values which actu-~ 


ally may occur. . 


The Flight-Phase Concept 


Similar arguments to those of the preceding section may be 
employed to justify ignoring the effect of the Earth and Mars grav- 
ity fields over certain portions of the flight. Certainly if a com- 
puter is designed to employ the values of these fields there will 
be a certain sensitivity level below which no cognizance of these 


values will be taken. 


Thus we may write with some degree of justification the 


following equations; 


[Roalt (eee Ca Phase II (3-14) 
(aval eae Says Phase I (3-15) 
[Ryaly = ft Gierk Phase III (3-16) 


where the limits of the Phases are to be identified. 


Thus by consideration of the trajectory and the magnitudes 
of the forces involved we reduce our problem from one of 12 bodies 


to three separate two-body problems. 
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TAS EE eal 


This table displays values for the parenthetic terms of 


equation 3-~¢ In the general form we have; 


a = Zoe ] 
G.,°¢ Gw)l= kh Sa = 
jA io j OJ Rp 2 R 2 
jA jo 
where k. is the surface gravity of the indicated body 


og is the radius of the body in AU 


Ria is the distance separating the body and the rocket in AU 


7 is the distance separating the body and point 0 in AU 


2.78 x 10° Pee Ome 


~32.2 —_—— ~32.2 


4 4 4 


lx 108 * 1x 107 * ) x10 


2 00Gb io=” | somos ame | Skuse: loe- 


7 fi it 


1.635 x10 6.21 x10. 6.2] x10 


~12.25 ~12.25 


fi HI 1 


7 lellom 9,21x10- 9,21x10. 


8 8 8 


4.025 x10 
10 


4.025 x10 
10 


12732 ca Om 


10 


4.65x10- Tt sole Tew ltrs Om 


10 


ZO <1 2 Ole 10. 





* Dy judicious choice of firing time 
t Assuming that Pluto and Earth have the same 
value for surface gravity. 


# Sample calcuiations are given in Table 3. 2. 
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Sample Calculation for the Values Listed in Table 3 1 
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CHAPTER 4 


THE TRAJECTORY 


A, SOME POSSIBLE TRAJECTORIES 


Basic Criterion for Selection 


Selection of the type of trajectory to be followed during the 
mid-course phase of an interplanetary exploratory trip of any 
type is dependent upon the basic criterion of when, in the future 
we contemplate attempting the flight. If we desire to make the 
attempt within less than ten years it is reasonable to postulate 
that presently developed power plants and guidance systems must 
be used with considerable refinement a necessity. If the attempt 
will not be made until some time after 1970, more exotic power 
plants and guidance can probably be postulated. However, the 
basic trajectories can be summarized and a selection made for 


further study in this paper. 


Classifications of Trajectories 


All trajectories, based on thrust available, can be divided 
into three classifications: (1) Powered, (2) Unpowered and (3) 
Combination of Powered and Unpowered. Although this statement 
is simple, the actual engineering and instrumentation re- 


quirements are exceedingly complex and somewhat dissimilar. 


Powered Trajectories 


The term ''Powered Trajectory" implies continuous thrust with 
either acceleration or deceleration used during the entire flight. 
Also implied is the ability to apply this thrust in any direction by 
multiple thrust units or vehicle reorientation with respect to the 


trajectory. 


The upper limit of this type trajectory can be considered to 
be strictly interceptor, while the lower limit is that of a distorted 
Keplerian Ellipse. Interceptor trajectories, using proportion- 
al navigation, collision course, line of sight, or constant Teed 
similar to present day antiaircraft and air-towir guided missiles, 
require nearly unlimited thrust. However, many simplifications 
in guidance methods can be realized due to reduced time of flight 
and reduction in accuracy requirements if high speed inter- 
ceptor methods are used. 


With small thrust power units such as the ion gun, solar 


furnace, or Nuclear-chemical, Big E196 


the trajectory can be 
a logarithmic spiral such as developed by Forbes or a series 
of constantly changing ellipses with the actual trajectory being 
the envelope oi the series. In addition the vehicle could be 
forced to traverse a predetermined elliptical path. This latter 


method would be extremely wasteful of both fuel and time. 


Powered trajectories are difficult to compute and almost 
infinite in number. The practicality of detailed powered 
trajectory study with any more than a general overall look 
at the field appears to be limited at presenti until the perfection of 


Suitable thrust devices warrant the effort. 


Unpowered Trajectories 


Basically the unpowered trajectory is a ballistic flight - 


Borat ork 
21, 44, 60 


icorporating the Fire Control Problem 
determining the path of a comet in space. This implies 
the application of a large impulsive thrust applied over a very 

Short interval of time with respect to total time of flight. This 
would result in a hyperbolic path with respect to a planet during 


a. e 
aes At the end of this time the hyperbolic 


application of thrust. 
trajectory would merge with an elliptical path with respect to the 
Sun, and the vehicle would be in free flight, acted upon basically 


by gravitational attraction of the sun with slight perturbations 
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6 


caused by the planets, (seeSection 3. 6). ; 

A trajectory of this sort leads to simplification or 
elimination of a guidance system for mid-course. However, 
extreme accuracy in the orientation of the velocity vector at 
termination of impulse is mandatory. Very small errors in 
magnitude or direction, or unforseen gravitational gradients 
during transit can cause large errors at the trajectory terminus 


as shown in Section, 8.6. 


Combination Trajectories 


Trajectories containing thrust periods and ballistic no- 
thrust periods are complicated by virtue of attempting to combine 
two Separately difficult problems into one. This method does 
permit slightly lower accuracies of the initial velocity vector 
but requires accurate mid-course position determination to be 


of any use. 


The basic advantage of this method is that presently developed 
thrust units can be used for correcting or changing the trajectory 
during mid-course, thereby increasing the probability that the 
mission will be completed. Using contemplated low thrust cor- 


ots ie powered phases could be of long duration and 


rective units 
be equally as effective as high thrust, short burning time chemical 


units. 


iMenyecvory oclecied tor Puriner onic, 


The basic mid-course trajectory chosen for further study in 
this paper is: that of a combination of unpowered and powered. 
The 1deal ov hoped for situation is one in which no corrective 
thrust is needed during the time of transit. This would mean that 
a relatively unperturbed ellipse could be flown within the limits of 
terminal accuracy requirements. However, if it becomes ap- 
parent during transit that corrective action is needed, a new 


ellipse is computed and proper thrust applied to get on the new 


45 


trajectory. The position determination accuracy requirements 
for such a system are rather high as are the thrust control 
requirements. However, it is felt that these requirements are 


not entirely beyond present day engineering capabilities. 


& 
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4.7 


B. THE UNPOWERED TRAJECTORY 


fatroduction 


To intelligently decide upon a trajectory to be traveled to Mars 
one must weigh many aspects of the problem. In previous por- 
tions of this paper we have investigated some properties which 
influence a vehicle in flight in the Solar System; and, in addition, 
investigated some vehicular behaviors and parameters which 


could be selected or adjusted to our benefit. 


As previously seantined the ''unpowered'' trajectory is a 
trajectory provided with an initial impulsive velocity, which, when 
perfectly instituted, needs no corrective action. This type tra~ 
jectory lends itself well to analysis since the vehicle is essentially 
a body acting in a central force field (subject to perturbations). 
Primarily due to considerations of available fuel and performance 
characteristics it was decided early in this work to proceed with 
an analysis of this type trajectory. It is emphasized here that 
this type trajectory may by no means be considered more at- 
tractive than any other because of any inherent property associ- 
ated with it. It simply appears more feasible due to the current 


state of the art in the development of rocket vehicles. 


ommiennics appears to be the first author (circa 1925) to 
seriously consider the solution of space flight problems. One of 
the most important results of his excellent study is the description 
of the minimum energy ellipse, which as the name implies re- 
quires the least rocket velocity of any configuration of two planets. 
This trajectory has become known as the ''Hohmann" ellipse and 
has the characteristic that the arrival of the rocket at Mars 


occurs when Mars is diametrically opposite Earth position at 


> 
cn 


4. 


8 


time of firing in their respective orbits about the Sun. Note 


Pae. 4. | 


ol, 92 agrees with Hohmann and shows some 


Preston-Thomas 
simple relations of the Hohmann ellipse with other ellipses. 
Other than this the literature shows a decided lack of information 


on other trajectories of the 'unpowered’ type. 


It is the purpose of this chapter to Show some relations exist- 
ing among the "unpowered' trajectories. These are the tra- 
jectories which when corrected for errors are described as the 


combination trajectories of Section 4.5 


Basic Relations 


The basic relations of a vehicle in motion under the influence 


44. 
of a central force field are set forth by Moulton -, Smart’? 
Kooy and Uytenbogaart?”, and others. Those general equations 


selected by the authors which appear to best fit the problem are: 


Z 
yo a(l-e ) (4-1) 


l-e cos 6 
L 
WY ad (le Stee (a2) 
tan @ we (4-3) 


Sl) Stereos so 


4. 
A fourth equation a which will be of interest is 


EN eee sin 9 ._-lle-cos@ m7 _ 2rt 4 
oa a a aera a = (4-4) 


l- e cos @ ue 
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A general trajectory situation is shown in plan view in Fig.4. 2 


Figs. 4.2 and, 3illustrate the quantities appearing in the 


foregoing equations. 


Hohmann Ellipse 





Fig. 4.1 The Hohmann Ellipse 


Some Simplified Trajectory Parameters 


We logically begin our analysis by considering the effect of 
the most massive member of the Solar System, the Sun. Re- 
ferring to Fig.4 2; for a given trajectory we may establish 
initial and final positional conditions (i.e. - we know we are on 


the Earth and hopefully we will arrive at Mars). Hence we write 


(4-1) as 


2 
_ al -e)_ 
"To 1-ecos 0, (4-5) 


and 


_ afl - 2). 


r, = (4-6) 
A l1-=ecos Oy 
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= 
rocket ellipse 


aphelion 


Fig. 4.2 General trajectory configuration 





Fig. 4.3 Some vector-angle relations for the general 
trajectory 
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where eee | ROG | at t= t,, and i -| Rod|at t=t) + te. 
Similar definitions apply to the angles 6. and 0 A’ 


Let r = 1 AU 
O 


luetr. = 1.52387 AU 


A 


Dividing (4-5) by (4-6) and rearranging, we have 


_ 0.5237 
° ~ 1.5237 cos, - cos @ (4-7) 
O 
Then from (4-1) 
(1 -ecos Dy ee 
OSE; in (4-8) 
l-e 
(4-9) 





where ais in AU. 


This defines two of the parameters (a and e) of any ellipse. 


We may now establish the vector velocity of the rocket wth 


respect to the Sun in terms of the trajectory parameters. 


Let Ne = V@ a at t= t- Then (4-2) may be written 


| V. = 97, 500 N Onene | ft/sec (4-10) 


where a is in AU. 


Let %, be the angle shown in Fig. 44 at t= t- Then 


(4-11) 








Vo@ 


TO SUN 


Pitas Velocities relative to Baris 





TO SUN 


Pig. 4.5 Velocities relative to Mars 


O00 
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Additionally, from Kepler's Laws 


Co A See (4-12 
T= NGM * 4-12) 


and the period of the ellipse is known. 


The boxed equations summarize pertinent parameters in the 
study of rocket motion with respect to the Sun, and form a basis 
for the calculations of velocity vector requirements with respect 


to Earth and Mars. 


It might be argued here that greater accuracy should have 
been demanded of these basic equations (i.e. - Use 1.523688 AU 
for | od | a Of 5230 Au ine a unit is 
known accurately only to about one part in 10° hence no good 
purpose could be served by extending the required accuracy of 
these equations beyond the fourth decimal, at least for a pre- 


liminary analysis. 


Calculations were carried out for about 50 elliptical tra- 
jectories and the values for e and| Vo] are plotted in Fig 4.6 
and 4.7 as aids in future studies. These plots presuppose that 
Ie and aN are given, known, or selected. A logical basis for 
the selection of O and On is the subject of another section. 
Eillipticity <1 was arbitrarily chosen since excessively high 
rocket velocities are required for parabolic or hyperbolic tra- 


jectories. 


Rocket Velocity Vectors with respect to Earth and Mars 


In the preceding section equations were developed to express 
parameters of the elliptical trajectory and initial velocity of the 
rocket with respect to the Sun. Obviously, since Earth and Mars 
are also bodies in motion with respect to the Sun we want the 
velocity of the rocket with respect to these bodies. The relations 
are Simple, geometric in nature, and are developed herewith. 

At t = ty 
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Fig. 4.7 Rocket velocity magnitude with respect to the Sun asa 


function of G and @ A 


as defined in Definition Summary 1 


From equation (4-3) 





e sin @ 
@ = “an 
O 1 - e cos @ 
O 
And from Fig. 4.4 
W 1 V_sin @ 
= tan Once ee eee 
Dir V., cos¢. - Ve 
O O 


V@a : COs Y@a 


ae VOM eet = 4 + 4, 


= e sin aN 





= to ———— 
Dai " We cos 6, 


and from Fig. 4.5 and equation (4-14). 


V, sin @¢ 
Von E tan™ = “7 cos @ 
= A A 





(4-13) 


(4-14) 


(4-15) 


(4-16) 


(4-17) 


(4-18) 


(4-19) 


(4-20) 





V, - V Ss 
O67 Vin.cO8A oat 


VOSA : cos Win 


If ry = 1.0000 AU andr, = 1.5237 AU then 


A 
Yo@ = 97,500 ft/sec 
and 
Vod = 78, 990 ft/sec 


Figures 4.8 and 4.9show Vey K and VIA as functions of a, 
and On: These values are important primarily as a step in 
solving the trajectory problem and are given here for possible 


future use by the reader. 


Hidmeror Flight 

The calculation of time of flight is necessary at this point in 
order to determine the angular separation of Earth and Mars at 
the instant of rocket firing. Furthermore, it is necessary to 
know t, in sidereal or calendar time since the inclination of Mars 


orbit with the ecliptic is a function of real time. 


If we substitute t. and t q 2s defined in Definition Summary 1 
each in turn in Equation (4-4) and recognize (= Q. and 


@ =0@, we have 
O d 


BN ae sin 9 ] e - cos @ - GMo 

Pane seas 0: ee | oe ae — : 

1 -ecos 9 ese 1 =-e cos @ T5 A 3 (4-22) 
O O 3 

and 
entmes sin @ e - cos 9 GM 
A _-l A T ea) 

l1=ecos @ Su 1 -ecos @ UE = Ye 3 (4-23) 

A A a 


Substituting te = ty - te after subtracting (4-23) from(4=22) we 


have 
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(4-24) 





where GM = 0.01716 if a is given in AU. t, is in days (sidereal). 


Values for a large number of trajectory time of flights were 
calculated and the results plottedin Fig. 4.10. Fig. 4.2 illus- 


braves pertinent factors in the foregoing equations. 





4.12 Angular Separation of Earth and Mars 
Having found the time or flight for a given trajectory we may 
compute the angillar separation of Hearth and Mars at the instant 
of firing (6. - @.). since Mars is in motion about the Sun(@. - 6.) 
is a quantity which must be known in order to fire at the correct 
time to produce an intercept. It is evident from Definition Summary l 


that 


a. =0, +6 st. (4-25) 


“7 
It may be shown” that 


: 5, (1 - e cos > s\e 
os) ae yp on 
= constant (1 - e cos Q$)" (4-27) 


Sec Go dsis a function of fogwe are inexorably engaged with 
a quantity which varies not only with real time but also with ang- 
ular position. A convenient method of expressing basty is to plot 
ores as a function of time; then, knowing te and choosing A we can 
obtain the angle 


mi = 6. - 0 (4-28) 
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and 
(6. - 6.) = (60 - Gy) - (6. - 64) (4-29) 


Tne aforementioned plot appears as Figs” 401) .. To assise 
in the logical selection of i an overlay (Fig. 4,12; in en- 
velope attached to back cover) has been constructed. The basis 
of construction and description of use is given in Information 


Summary 4,1 
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MARS ANGLE IN DEGREES FROM EARTH APHELION 


JUTE 
JA 


240 
Mars 
Aphelion 


Mars 
Perihelion 
. - 


0 O02 04 06 O08 $%I10 jJ2 j%4 16 18 20 2.2 
TIME IN EARTH SIDEREAL YEARS 


Fig. 4.11 Mars angular position measured from Earth aphelion 
as a function of time 
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Fig. 4.12 (Overlay #1) 
Earth angular position measured from earth 


aphelion as a function of time 


(In pocket on back cover) 
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INFORMATION SUMMARY 4. 1 


Construction and Use of Overlay #1. 


Assuming the "state of the art" of interplanetary travel is 
such that an Earth-Mars trip is planned to occur about 1980, 
we note’ i that on January 1, 1979 Earth is 177° ahead of her 
aphelion (see insert on Fig. 4.12 ). Mars is 2° behind 
Farth's aphelion. In other words Mars is 181° ahead of the 
Farth. During the succeeding 1.135 sidereal years Earth 
overtakes Mars, and*they are in conjunction DOT. past Barth 
aphelion (about February 17, 1980). Hence from January 1, 1979 
to February 17, 1980 oe = 6.) takes on all values from 0 to 7. 
With all pertinent parameters of both Earth and Mars orbits 
ne wie Fig. 4.jiand 4.12 were drawn so that information for 
the period in question would be easily accessible. The as- 
Sumption has been implicitly made herein that Mars aphelion 
precesses only slightly about the Sun (a few seconds of arc in 
about 10, 000 years). oF 


(a) Basis for Construction of Overlay 


The angular position of Earth is indicated for January 1,1955 
and will repeat every sidereal year. Mars angular position is 
marked on graph 4.11 for January 1, 1955, but its period is not 
an even multiple of years, hence its angular position for any 
given time must be found either from an Ephemeris or by the 


following method: 


(Page 1 of 3) 





Ge 





Divide the difference in time between January 1, 1955 and 
the time for which Mars angular position is desired by the 
number 1.8& (Mars period in sidereal years). Slide the mark 
for January 1, 1955 on the overlay to the right of the mark for 
January 1, 1955 on the underlying graph an amount equal in 
sidereal years to the remainder in the foregoing division. The 
real time scales have now been superimposed so that the 
angular positions of Earth and Mars are coincident on the time 
scale. The values of or and GO. may be picked off the curves 
and ordinate scale. If a time other than January 1, for a given 
year is desired, make auxilliary pencil marks for Earth position 


and use that rather than January 1, as a reference. 


(b) To find (@. - 6.) for a given t. 





During the period of interest all values of (0, - 6.) from 
Oto m occur, hence set the overlay as before and read te corre- 
sponding to the given ss and O. as a decimal fraction of a year 


after the chosen reference time. 


(c) To find t. for a given @_ and @ 
O O S 


This problem is somewhat more complex. We are primar- 
ily COncerned with a real time firing time in late 1979. TEhis 
is true because it is about this time that (6. - 6.) takes on 
values in which we are interested (note Section 4.12 and Fig. 
4.11 and 4,12 ). Fortuitously; Mars' angular rate in late 
1979 is approximately constant at 0.47 degrees per sidereal 
day. If we multiply 0.47 times the te this determines a rough 
value of (6. - Oy): Subtracting (6. - Gy) from (6. - GA) pro- 
vides (6. -- 6.) for which a t, can be found as explained in part 
(c). By nothing the angular positions of Earth and Mars and by 
moving a straight-edge up or down the superimposed curves 
one may obtain a more accurate value of (@- 6.) for the given 


data. 
(Page 2 of 3) 


= 
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(d) To find t, for a givent,, @ and 0, 


January 1, 1955 was chosen as reference time because 


accurate values were available to the authors. 


The method described herein is a rough graphical method 
of finding required quantities for preliminary work. Actual 


calculation would of necessity be made by computer. 


(e) Notes 





(Page 3 of 3) 





4.13 Ven and Vo A as functions of (0. = 6.) 


Since we now may find (6. - 6.) for any given trajectory, 
it appears pertinent to note any relations which may exist 
between this angle and the velocity required to fly the trajecto- 
ry. From an inspection of the basic equations in Section 4.8 
it is obvious that the relationship is quite complex. For this 
reason Von and GUN were plotted graphically for a large 
number of trajectories in a manner analogous to that of 
Preston-Thomas*-. The minimum Vo a was found to be the 
envelope of the curves as 6, took on various values. (See Fig. 

4.13 ) The minimum Vd a Was found to be identically the 
curve for On - 0°. (See Fis. 4.14 ) This is not inconsistent 
Since the first case represented minimum velocities for a 
non-maximum radius vector of an ellipse; whereas the second 
case represented the minimum velocities that exist at aphelion 
(maximum radius) of an ellipse. In both cases the Hohmann 
Ellipse (6, = 180°; 6,= 0°) provided the absolute minimum 


velocities. 


4.14 The Rocket Plane (A plane) Concept 


Let us define the plane containing Rog and Waya at t= vs 
as the rocket plane (A plane). If the rocket were not perturbed 


during its flight it could remain in this A plane (this being a 


op) 
Cn 


(Vea)t=t, IN THOUSANDS OF FT/SEC 
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Fig. 4.13 Rocket VEU hay magnitude with respect to Earth as a function 
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property of the motion of a body ina central force relearn 
The A plane need not necessarily, anc generally will not,lie in 
the Earth ecliptic plane since Mars passes through the Earth 
ecliptic plane only twice every 1.86 sidereal years. Unfortu- 
nately, for calculation purposes, the orientation of Von in 
Solar System space depends on the real time firing time and 
time of flight. Hence every possible trajectory requires a 
complete individual solution. The labor involved is so moun- 
tainous that solution by high-speed digital computer is mande- 
tory; and even these solutions must be based ona selected 


Mitre wane time, 


If we can show that the displacement of the rocket normal 
to the A plane is negligibly small we may be justified in saying 
that for practical pusposes the rocket remains essentially in 
the rocket plane and our problem reduces to two dimensions. 
Neglecting the possible perturbations of bodies such as meteors, 
distant planets, asteroids, etc., the bodies causing displacement 
of the rocket from the A plane are Earth, Moon, and Mars. In 
an effort to show the minor nature of the aforementioned dis- 
placements calculations were made ona single ''typical'' 
trajectory. This trajectory was chosen because oe =10, 000 
ft/sec, which is a velocity probably within the capabilities 
achievable in the near future Since we are concerned with 
the period from late 1979 to early 19360 the firing time is es- 
tablished by Vor ~ oe) to be about December 6, 1979 Note 
Passe i> and 4.12. " 


Referring to Derivation Summary 4.1 the cumulative dis-~ 
placement of the rocket from the A plane was found by ac- 
cumulating the values determined from Equation 1] over the 
trajectory traveled by an unperturbed rocket(see Fig 4.19). 
Certain errors are inherent in the calculations but will be seen 


to be negligibly small. The arc length of the trajectory 


ROCKET ELLIPSE DATE OF FIRING 
APHELION APPROX. DEC.6,1979 
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Fig, 4.15 Angular quantities of a representative 
trajectory discussed in section 4.14 
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Fig. 4.16 A typical trajectory showing the effect of tilt to account for 
non-co-planarity of Earth and Mars motion 
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AN IMAGINARY CYLINDER 
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Fig. 4.17 Assisting supplementary view to figure 4. 16 
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DERIVATION OF THE EQUATION FOR THE DETERMINATION 
OF THE DISPLACEMENT OF THE ROCKET NORMAL TO THE 
“ROCKET PLANE" a: 


Four relations are available from consideration of the geo- 
metrical properties of the trajectories or orbits. Two assumptions 
are made which greatly simplify the solution. The error introduced 
by these assumptions is on the "high" side and will be shown to be 


negligible (Fig. 4.16 refers). 


eG gsin & (1) 

Wee (i ax sin (A + 6, - Te) ) (2) 

ZA = Ro, sini, (3) 

4 rcs sin id (4) 
From (3) and (4) 

4 © ain = sin i (5) 

A Ro, o) 


(a) Basic Equations 
To aid in clarification of the development that follows and 
the assumptions made, the following list of symbols is presented: 


(Sorthj Magnitude of that component of G, normal to 


the "rocket plane”. 


DERIVATION SUMMARY 4.1 
(Page 1 of 4) 


1 Inclination of indicated quantity with reference 


J 

plane, or with ecliptic plane if no reference 
is given-positive direction in the conventional 
sense. 

(i ). Maximum inclination as defined above. 

max’ j 

Z. Displacement of indicated quantity normal to 
indicated reference or to ecliptic plane if no 
reference is given. 

A Angle as shown in Fig. 4.16 

& Eccentric anomaly as shown in Fig. 4.16 

5 Increment (Magnitude as defined). 


(b) Definitions 


Referring to Fig.4.16 and4.17the two assumptions made are 
(1) that the "rocket ellipse'' is almost circular about its 


center (C). 


(2) that the line CC' in Fig.4.16 is negligibly small. 


(c) Assumptions 


To derive the required equation we solve ( 3 ), then 
( 4  ) for the final condition a5= ae We then find iA from 
Geo.) and eee) from (1 ). Having discovered eed 


for a given trajectory we are now able to write iA as a function 
of ve) (Equation 1 ). 


Any [ z may be found thusly: 


Aplane 


DERIVATION SUMMARY 4.1 
(Page 2 of 4) 





~) 
Cj! 





Poo ric. 4.17 ; 





4 pplane! ats Bron oe Yl at 
and : 
ZS = Rye) sin Slay 
obviously 
[4d] Aplane 2 'Nolane! ats = 2G 
now 
ee [2d] Aplane 
wa 
PS | 
measured from 
rocket 
then 
ho eee Ze) A 
. ~_ * Oplane 
l Sorthld - JGeyatiey. rp) PF ral IN 
and 


[ 6Z (due to Mars) 


A Aplane 


orthld 2 
Ey a (6 t) 


y 


2 5 
. l Sorthle Co T)mex 7 lon 110 


if 6t is in sidereal days 
Eorth is in ft/sec? 
6 ZA 18 in miles 


DERIVATION SUMMARY 4.1 
(Page 3 of 4) 


ne er teen Tig ae ee 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 








A similar derivation may be made for the cis-lunar 


Barycenter. 


(d) Derivation 


DERIVATION SUMMARY 4.1 
(Page 4 of 4) 


~.] 
=~] 


traveled by the rocket is approximately 354, 000, 000 miles; the 
cumulative displacement from the A plane is approximately 5, 500 
miles: a ratio of about 1.55 to 100,000, certainly negligible. 
Further, the calculations show that the displacement due to the 
mass attraction of Mars is only 33 miles, the remainder being 
the effect of the Earth-Moon mass attraction. About 3, 250 

miles displacement occurs in the first five days of travel; hence, 
if a computer program were established to counteract this, the 


displacement could be reduced even further. 


The obvious advantages of reducing the problem to two 
dimensions are many-fold. These include (a) reduction in 
weight and complexity of the guidance mechanism; (b) simpli- 
fication of component construction; (c) savings in fuel; (d) im- 
provement of payload ratio; (e) saving in material, labor, and 


cost; (f) contraction of production time-tables; and many others. 
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CHAPTER 0 


DEPARTURE AND ARRIVAL 
ORBITS 


Introduction 


The capability of a rocket to achieve required velocities de- 
pends so critically upon the energy carried within itself that 
every ellort must be made 10 minimize these energy requiremenrc: 
This section investigates one aspect of this minimization -- the 
selection of departure and arrival orbits about the respective 
planets to minimize kinetic energy, hence velocity, that the 
rocket must provide itself intorder 10 fly 4 Sivem trajecciom: 
lees OF | illustrates the relationships of pertinent velocities. 


For complete definitions see Definition Summary 1 


am 
a 
Pe a, Ce ee eee = 
\ Te one 
\ — “pa @A)g 
\ NT inc '@s 
p Ney = 
Og = _ . = Vi@A)s 
saa ok V % _ 
<4 OBA V 
V @s 
BA i \ > 
y 
Yoo © 


Fig. 5,1 Pertinent Velocity Relations 


Seeeeoptimum Orbit Radi 


In order to fly any given trajectory, or ellipse, a velocity of 
the rocket with respect to the Sun (Vo A) must somehow be attained. 
Referring to Section 4.9, we note that this is the velocity vector 
required for an intercept of the orbit of Mars, treating both Earth 
and Mars initially as massless points. However, this is obviously 
not the case. Referring to Fig. 5.2 we can say, figuratively, 
that we must climb up the side of the Earth potential "well", over 
the intervening "hump', then down into Mars potential "well. 
Jupiter's potential ''well" is shown for comparison. From Fig.5.1 
we note that the rocket already possesses a certain velocity with 
respect to the Sun by virtue of the Earth's orbital motion, thus the 
rocket must be provided the velocity, Ve a: lmorder to "over- 
come’ the effect of the earth gravity field, the rocket must be 
provided a eyodiee velocity, Vi@a)e’ The rocket also possesses 
a velocity, VQ ae by virtue of being in orbit about the earth, 
therefore it remains for the rocket to provide by its own efforts 
an increment of velocity Va Initally Von will be along the 
dotted line labeled in Fig. 5.1, however, the earth gravity effect 
decreases rapidly so that VO A HN almost identically on the solid 
line. In addition the ratio "Te 5/Ro@ is of the order 10 “, hence 
the displacement of the rocket from the earth 'O> is negligible 


for purposes of calculating initial velocity requirements. 


Regardless of validity, if VOa is known we may obtain 
relations between a and T@ 5° This has been done in Derivation 
Summary 5.1, Although the relations have been developed for the 
planet Earth the derivation is perfectly general and applies to any 


body; with different constants. 


Parts (b) and (c) of Derivation Summary 5.1 develop and 
prove the optimal relationship of V,_, and TQ: Part (d) expresses 


the relation non-dimensionally. 


Pertinent ranges of Vi vs t@ « have been plotted in Fig. 


c 
9.3 and the non-dimensional function is shown over an applicable 


mange in Big, 9.4. 
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DISTANCE IN MILES 


Fig. 5.2 Velocity of "escape" from the Sun and some planets 


as a function of distance 


10 


60 


w aon 
So So 


a 
So 


=o 
(Tislopy IN MILES x 10 


20 





(Vine)opt IN FT/SEC x10°° 


Fig. 5.3 Optimal relationships between Mec and ris for 
Earth and Mars 
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is the frequency of occurence of the rather special Earth-Mars 
configurations required bythe Hohmann ellipse. This, among 


other considerations has been taken up in preceding sections. 


4 


= 


Use of Optimum Orbit Concept 


Almost paradoxially, the optimum orbit concept is of 
particular use in the selection of non-optimum orbits or tra- 
jectories. As an example, assume that a rocket has been 
developed with a capacity to provide itself witha err of 
approximately 10,000 ft/sec. The optimal orbit radius corre- 
sponding tc 10,000 ft/sec, as an optimal velocity, from graph 5 3 
iS approximately 27,500 miles above Earth center. If 27,500 
miles is unacceptable due to other considerations we can reduce 
(from graph 9.4 ) this radius by 42.5% simply by allowing a 2% 
increase in the required velocity. In other words, for Vos 
iG, 200 ib)/Ssec. 7: @s > 15,800 miles, or approximately 12, 000 
miles above the Earth's surface. It must be borne in mind, 
however, that these relations exist only for that trajectory for 


wineh (V. -) .. = lOsG00 1h see 
ine opt 


Applying the foregoing reasoning to the 'Hohmann'" type 


ellipse -(V.__) = 6, 818 ft/sec; (r ) = 57,300 miles - we 
inc’ opt Bs‘ opt 

may reduce lg, to 33,000 miles by increasing ao to 6, 956 

ft/sec. Hence. if we can provide the energy to increase Velo 

by 18 ft/sec we not only gain the obvious advantages of a lesser 


Y o> but also the flexibility inherent ina greater freedom of 
choice in the angular positions of Earth and Mars. This has been 


discussed more fully in Chapter 4 


The Least Orbit Concept 


If the problem is posed that V,_, = 7, 500 ft/sec. , what is 
the least orbit the rocket may depart from and still travel a 
trajectory to arrive at Mars? This obviously is not an optimum 
orbit problem as previously posed. For (V.__) = 7,500 ft/sec 


inc’ opt 


(r alloy oh = 47,000 miles (graph 5.3). For Ee ont = 6, 818 


ft/sec (the minimum 'Hohmann" ellipse velocity) the ratio Vee 


Dene ope 1.1 and from graph 5.4 the ratio r Bares ®s)opt 


= 0.327 andr Qs" 18,750 miles The limiting restriction herein 
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DEPARTURE AND ARRIVAL 
ORBITS 


Optimum Planetary Orbits 


In order to' 


‘escape’ any bodies mass attraction a vehicle 
must attain kenetic energy at least equal to the potential energy 


inherent in the field or fields to be "overcome"; hence 


KE = my | 1G; | dR (1A) 
or 
o GM 
SE =) dR (1B) 
OID 58 ee 
Following the convention @ = 1Q | we integrate 
GM. 
ee ian 8 (2) 
A jA 
but AE = 1/27. (3) 
A J 
2GM. 
hence Wee = ra (4) 


It has been Shown” that a rocket in circular orbit about a 


body possesses a so-called circular velocity. 


| DERIVATION SUMMARY 5.1 


(Page 1 of 5) 
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; V. 
Je Js 


From Fig. inthe case of the earth and a rocket 


Vi@Aje = V : + V, 


(7) 


Purther, by equating tne kinetic enersy requinemenrcmion 


the rocket we have. 


aN 


2 an 25 isp V eon, Wie 
@e 


Ap 4 


A® (DA)g 


(8) 


where the first term represents the kinetic energy requirement 


ponachieve the velocing V@ a? the second term preresents the 


kinetic energy requirement to escape’ the earth gravity field; 


and the third term describes the sum of the first and second. 


From (8) we have 
ae 2 
Hence, substituting in (7) and rearranging 


. 2 


shale: ISS 
and if initially R ee then 
Rae ~ 2GM, [G GMep 
Ve Te ee 
inc SP) ea De 
Ds S 


(a) Development of Vinc 12 terms of rocket orbital radius 


DERIVATION SUMMARY 5.1 
(Page 2 of, 


of” 


(9) 








In order that rg@. be the optimum radius on a kinetic energy 
basis V,.., must be a minimum for 1 (TQ s)ont which is the 
equivalent of the conditions . 


ic 0 
(1) ae 5 = 4() (12) 


brg c: 


To satisfy (1) we take the first partial derivative of (11) 
and equate to zero. 
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To satisfy condition (2) above we take the second partial 


derivative of Masi with respect to 1D ¢ and evaluate at '@q 


c 
= @>5) opt to show it is indeed positive. 
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which is positive for positive values of the radical 
(c) Proof of Minimality of Qs) opt 
Dividing (11) by (7B) and noting (17A) we have 
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CHAP Tih. 


THRUST AND ENERGY CONSIDERATIONS 


A. THRUST REQUIREMENTS 


Introduction 


It is of interest to examine the thrust requirements and 
payload ratios for a rocket departing circular orbit about earth. 
such figures appear to be more or less well known for the 
journey from earth's surface into orbit (as for earth satellites). 
However, the orbit-to-orbit requirements are not commonly 
available, although easily calculable. This information is con- 
sistent with the optimal orbit concept developed previously and 
verifies the validity of such an approach for the conditions 


assumed. 


Although the data have been based upon the departure and 
arrival orbit concept, an extension to surface departure can 
easily be made by extrapolation to lower altitudes. The penalty 
of departing a lower altitude will be evident after an investigation 


of this section 


The Optimal Payload Ratio 


The change in velocity necessary to affect transfer from 
orbital motion about earth to departure path has been designated, 
Were. This change in velocity is accomplished by applying thrust 
for a controlled period of time. [f the power-on period is suf- 
ficiently brief. the maneuver trajectory can be neglected and the 
transition from circular orbit to departure path thought to occur 
inStantaneously at tangency. The justification of this asSumption 


will be discussed shortly after the development of the equations. 


The commonly used expressions for an idealized single 





step thrust study of this nature are derived from the equations of 


motion written along the orbit flight path. ee 


In circular orbit, however, for the above assumptions, D 


is presumed negligible and 6 = 0. 


a 
Mee 2 dt 


where M, = initial gross mass 
AM = mass of propellant 


t, = time of burning 


b 


A constant thrust device is assumed of specific impulse, I. 





M. 
b 
i © Me 5e 
dv _ _b D where =“!2 
dt M -M 1-é6t M, 
Py, t 
th b 


Integrating between t = 0 andt = th 


a Ni ea et (1- 5)= - Ig ln (PR) 
M, -M 
PR = Payload Ratio = a 
O 


It is to be remembered that this, V. applies only to 


INCY® ” 
the departure maneuver, as does PR@ - FOr Gilby miirore! Gellar 
orbit about Mars, a similar expression is obtained. The total 
payload ratio( PR,) for an orbit-to-orbit maneuver, ignoring 


corrective thrust requirements, is 


<a) 
itm 


PR, = PR» PR 


These equations have been expressed in tabular and 
graphical form for various values of specific impulse, which 
has been chosen as the propellant parameter. The subscripts 
in the headings of Table 6.1 refer to the three assumed values 
of specific impulse. Table6.1lforexample, lists Payload ratios 
for various sample trajectories developed by the methods of 
chapter 4 . The importance of improved propellant performance 


is noted. Although the corrective thrust requirements have not 


PRe 300 _*@350_ 
063 | 


Oi .139 


been included, reference 47 
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TABLE 6.1 
SOME REPRESENTATIVE TRAJECTORIES WITH CORRESPONDING 
PAYLOAD REQUIREMENTS 
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has indicated that a suitable payload ratio for corrective thrust 
requirements is .79. If the above figures are altered by this 
multiplier, it is seen that the most optimistic Payload Ratio to 

be hoped for is .ll. The calculations of Table 6.1 presume de- 
parture and arrival orbits at "optimal’' altitude as defined 
earlier. Since these altitudes inevitably vary for various 
trajectories, little correlation is to be expected between Payload 
Ratio (PR) and angular separation at instant of fire (6. - 6.) 
other than to indicate a lower bound for certain specified para- 


meters, 


6.3 An Interpretation of Optimum" Altitude 


At first glance, it appears that the higher the optimal 
departure and arrival orbit altitudes, the more favorable the 
total Payload Ratio. To invalidate this idea as well as to demon- 
strate the advisability of employing optimal altitude two tra- 
jJectories have been selected for study: one is the familiar Hohmann 
semi-ellipse trajectory: the other merely a representative one 


selected for contrast. Both are defined by the * in Table I. 


Ifthe PR@ and PR, for these two trajectories are 
calculated as functions of altitude (including optimal altitude) on 
a non-dimensional basis, as in Fig.6.land6.2 it can be seen 
that a peak paylcad ratio is achieved at precisely(r OPT’ A slight 


fall-off occurs at values oft >1. Thus, the choice of 
° OPT 
optimum altitudes, in the sense defined, insures a maximum 
total payload ratio for any given specification of Vine’ and 
Vine! . The great saving in Payload Ratio is noted when one 
=> 


r 
departs from circular orbit altitude for ee 2.0. 
s°OPT 


6.4 The Assumption of Impulsive Thrust 





A justification of the assumption of impulsive thrust is 
dependent upon a comparison of ern time and period of orbit. 


Obviously, for a low thrust device of considerable duration, the 





Cis 


(Tis )opt 


Payload ratios VS. orbit altitude for the Hohmann ellipse 


6.1 


Fig. 
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Fig. 6.2 Payload ratios vs. orbit altitude for a representative trajectory 
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Fig.6.3 Orbit period as a function of altitude above planet center 


assumption is not valid. 


The period of orbit for both Mars and Earth is plotted in 
Fig.6.3 as a function of distance from Earth and Mars center. 
If we assume that instrumentation is capable of withstanding an 
axial load of 10g upon departure from Earth ofbit the minimum 
acceptable burning time can be determined. oe Careful design 
can then be employed to produce a favorable burning time. The 
minimum burning times as functions of orbital period for optimum 
Earth and Mars orbit are listed in Table 6.2 for the trajectories 


previously considered: A specific impulse of 350 seconds is 


assumed. 


T (min) 
"S ee "so ie 





TABLE 6.2 
A COMPARISON OF ORBITAL PERIOD AND BURNING TIME 


Sie 
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B. ENERGY CONSIDERATIONS 


Introduction 


It would seem appropriate at this stage to mention opti- 
mization of the trajectory. The term, optimization, however 
invariably leads to ambiguity and requires explanation. Do we 
desire optimization with regard to time, fuel expenditure, or 
initial velocity? -- or possibly cost? Previously, an attempt was 
made to minimize the initial velocity requirements, but it is 
clear that this does not necessarily represent the optimum as 


far as other variables are concerned. 


An overriding consideration during early interplanetary 
flights will undoubtedly be minimization of fuel expenditure. By 
minimization of fuel expenditure, we imply a determination of 
not only the minimum duration of powered flight but also the 
direction in which this power or thrust is applied. Accordingly, 
these are all considerations in the minimization problem. It is 
the intent of this section to describe the problem in general terms 
and summarize, without proof, the findings of other authors who 


have written at length on this particular subject. es Oo 


General Theory 


Let us first consider the problem of optimum transfer (with 
respect to fuel expenditure) between two terminals lying in a 
plane. If we consider an inertial orthogonal coordinate system, 


the equations of motion for the rocket are 


oe G,, = -cOS a 77 oe (6-1) 
YG - cose aM 6-2) 


A + G, = -cos 1a a (6-3) 
where G,. Gy, and ee are the components of ul resultant at- 
tractive force per unit mass and a, Band y are the angles 
between rocket thrust direction and coordinate axes. The ex- 
pression aa o is the well known simplified term for rocket 
thrust and is defined in the section describing thrust consider- 
ations. c is the effective exhaust velocity and M is the mass of 


the rocket at any time t. 


G, Gy G, depend explicitly on x, y, z, andt only. 
Squaring and adding 6-], 6-2, and 6-3. 


ie cee 7 a 2 ce Collie 2 2 Z 
[X+G,] +[ ¥+G.] +[Z+G,] =([yi Gi 1 [cos et+cos B+cos y] 
(6-4) 
Integrating, 
M ‘] ; 
O ee ye 2h ae Zi oye 
@ eM, = \ l x+ G,.] 15 aie Gy] sal erect G,] dt 


It is assumed that the boundary or terminal conditions are 


specified 
e > e Cf 
<— (t.) = X, x (to) = X, X, (t,) = Xx, x (t,) = Xx, 
e # : e 6 
oe (t.) = ee Y (to) = oe Yj (ty) = ag) Y) (t)) = Y 
® ® @ 6 
Zo i) = Zo Te (to) = 4, Z,(t,) = 2) ci 2) (t)) = aa 


The objective is to select the functions x(t), y(t), and z(t) so 
that the fuel EURO is a minimum (i.e., we desire initial fuel 
mass, M, as small! as possible and still attain our objective 


subject to the boundary conditions). 


The minimization of the above integral is a problem in the 
calculus of variations. Its solution is highly detailed and is not 


presented here. The conclusions and interpretations of existing 
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solutions are helpful and interesting, however. and are stated 


here without Boor i 


These conclusions are due to D. F. Lawden, one of the 


more active writers in the field. 


bet us consider a rocket orbiting the eantneci ne mean hy eee 
optimal path is desired for travel from this orbit to entry into 
circular orbit about another planet. The conclusions of Lawden 
indicate that an optimal path can only be followed in a vacuum by 
the application of certain impulsive thrusts from the rocket motor, 
the rocket coasting freely under gravity during the intervals 
between Such thrusts. These coast periods are essentially null- 
thrust arcs and are mihimal in the sense that neighboring arcs 
can be followed only by an actual expenditure of fuel. To establish 
an absolute minimal mode of transfer for a definite planet con- 
figuration, cases for various numbers of impulses must be com- 
puted’ and the minimum selected. Lawden suggests that two 
impulses usually lead to the optimal solution. It should be re- 
membered, however, that there will be a large number of optimal 
solutions corresponding to different planet positions. Only one of 
these will be an absolute optimal solution for all planet con- 
figurations. The absolute optimal solution for an Earth-Mars 
flight has been known theoretically for some time if certain 
fundamental assumptions are made. The analysis“~ is that of 
Hohmann, who concluded in 1925 that the problem of transferring 
a rocket from one circular orbit into another coplanar with the 
first and about the same attracting body would be most economical 
if effected along an elliptical trajectory tangential at its apses to 
both circular orbits. Fuel would be expended rapidly at each apse 
to carry the rocket into and out of the transfer ellipse. A state of 
free flight would exist between the two impulse trajectory. This 
classic development, the so-called Hohmann ellipse, has been 
verified by Lawden am is consistent with the remarks above. Note 


ihai the optimum path as developed and computed in this paper on 
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the basis of a minimum initial velocity is in fact the Hohmann 
ellipse. The importance, however, of knowing the optimum for 
other positions of Mars and Earth so that flight can be initiated 
often, is not to be underestimated. 


® 
Mars 


Orbit Rocket Path 





B A 
Decelerating Accelerating 
Impulse Impulse 


Fig. 6.4 The Theoretical Mode of Transfer 


6.%7 Rocket Thrust Direction 


Ordinarily, the mode of directing the thrust which results 
in maximum total energy is found to be alignment of the effective 
exhaust velocity with the velocity vector of the rocket (i.e. align 
the direction of thrust with the tangent to the trajectory. This 
will ensure that the rate of increase of total energy is maintained 
at its maximum value at every instant. This program will be an 
economical one since the object of the maneuver is to raise the 
total energy as rapidly as possible in order to achieve escape 
and sufficient velocity. However, it has been gaevedl that this 
program may not be an absolute optimal method, although close 
to it. Under some circumstances, when departing a circular orbit, 


it is not economical to maintain the rate of increase of energy as 


1G2 


its highest value over the entire trajectory of escape. The 
author claims that after having escaped, it is advantageous first 
to direct the thrust so as to oppose the motion. The rocket then 
drifts toward the center of attraction along an elliptical arc. At 
its closest point of approach, an impulsive thrust is generated in 
the direction of motion. The velocity is now large, and the total 
energy increases at a rapid rate. The rocket is transferred into 
a hyperbolic orbit along which it coasts to its destination and the 
overall fuel consumption is found to be less than if escape had 
been achieved by application of a single tangential thrust. The 
Beaders relerred 10 reference” for information as to how 
advantage can be taken of such maneuvers for general conservation 
of energy. A good approximation to the optimal programming of 
thrust direction under departure from circular orbit may be 
achieved by maintaining alignment between directions of thrust 


and. motion. 


Other writers i consider it an empirical rule that a journey 
from one to another stable orbit round a central body will be 
completed in the shortest time only if power is applied in two 


impulses, one at beginning and one at end of journey. 
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CHAPTER 7 


NAVIGATION IN SPACE 


A, POSITION IN SPACE 


General Discussion of Problem 


Mid-course interplanetary navigation differs quite markedly 
from terrestrial position determination due primarily to the need 
to obtain a fix within a volume instead of upon the surface of a 
sphere (Earth) of known radius. Terrestrial positions can be 
obtained by simple spherical trigonometric relations and can be 


further simplified to a problem in plane geometry and time. s 


Although it is obvious that the problem of Position in Space 
is three dimensional, several writers have postulated that a suf- 
ficiently accurate solution can be obtained by considering the 
planets to be coplanar and working with only two dimensions in 


the plane of the ecliptic. SP he 


Part of the argument for this 
approach is found in star tracker limitations, errors in astro- 
nomical measurements, the assumption that terminal guidance 
will correct for errors incurred during mid-course transit, and 
the relatively small angles of inclinations of the planet orbits 


from the ecliptic plane. 


Factors Indicating Necessity for ''On Board'' Automatic Solution 


if we postulate that accurate positional data is needed during 

mid-course transit, several factors indicated that the problem 
must be solved on board the vehicle and by automatic computers 
of some sort. 

1. Due to high vehicle velocities (in the neighborhood of 

1,200 miles per min) time delays of any sort should be 

avoided wherever possible. This eliminates manual 

computation and sight taking in the vehicle itself or on 


the earth. 


105 


2. Radio frequency transmission in space has a rather 
sizeable time delay. (A 10° mile round trip signal requires 


approximately 18 minutes) 


3, Radio frequency signal interference in space is un- 
certain and apparently unpredictable ‘when sent to or 
from Earth due to the Earth's atmosphere and extraneous 


radio frequency waves in space oot 


4. Power requirements for dependable radio trans- 
missions during mid-course appear to be far too great for 


installation on any vehicle within the near future. 


5. Optical tracking of the vehicle from earth is not possi- 


ble with any foreseeable astro telescopic equine a 


6. Radar tracking of the vehicle from earth is ruled out 
by 2 and 3. 


7.3 Selection of a Reference Coordinate System 


Celestial bodies may be located with respect to several 
different reference systems which are shown inFig. 7.1 and 
Table 7.1 . For mid-course interplanetary position de- 
termination two of these systems can be discarded immediately. 
These are the Geocentric Horizon and Equatorial. Coordinate 
systems using Galactic quantities are not particularly convenient 
for travel limited to our solar system. Of the remaining only the 
heliocentric ecliptic system will be considered for use in this 
paper. This system has the following features which lead to 


this choice. 


]. The reference zero latitude plane includes the Earth 


and Sun at all times. 
2. All other planet's orbits lie close to this plane. 


3. NumerousEphemerides are available based in this 
system simplifying programming of any memory drum for 


a computer 
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Establishing the Reference Coordinate System 


It appears that the reference system orientation will have to 
be obtained by use of the fixed stars, an extremely sensitive gyro 
inertial platform, i or a combination of the two. The use of 
gyroscopic elements by themselves for establishing and main- 
taining a chosen reference coordinate system to a sufficient degree 
of accuracy to determine accurate angular measurements and rates 
in space appears to be beyond engineering capabilities at present. 
All instrument gyros have drift uncertainties to some extent and 
very definite lower limits of sensitivities. To be useful for ac- 
curate interplanetary guidance the threshold of sensitivity would 
have to be in the order of 1078 9. Since it is necessary to con- 
template transit times of over 100 days, it appears that star 
trackers locked on the fixed stars will have to be used to a cer- 
tain degree in either assisting a gyro system or replacing it 
entirely. At certain periods of the transit it may be advantageous 
to use gyroscopic elements as a memory circuit or coast device 
while a star is temporarily obscured or its line of sight becomes 
nearly coincident to that of a stronger source of light such as a 
planet or the sun. It also appears that during periods of radical 
maneuvering or high thrust a gyroscopic device may be needed 


to insure holding track on a desired star. 


If we postulate the use of photoelectric star trackers, es- 
tablishing the desired reference coordinate system can be ac- 
complished by using only two fixed stars approximately 90° apart 
in the same meridian. However, three stars near quadrature to 
each other, or three stars near the chosen zero latitude spaced 
approximately 120° apart, will establish a more stable and secure 
reference. Monitoring three fixed stars would allow for the 
possibility of one of them becoming momentarily obscured or 
temporarily loosing tracking signal. Section 8,64 discusses 


the accuracy requirements for such systems. 
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Determining Celestial Latitude and Longitude of the Vehicle 


After establishing the reference coordinate plane and meridian, 
the celestial latitude and longitude of the vehicle can be determined 
from the line of sight to the sun as seen from the vehicle. This 
is illustrated in Fig. 7,2 and as shown, the vehicle heliocentric 
latitude is equal to the negative of the angle between the reference 
zero latitude plane and the vehicle sun line measured in a meridian 
of the reference coordinate system. The heliocentric longitude of 
the vehicle is equal to the angle between the reference zero longitude 
and the meridian containing the vehicle sun line measured in the 


direction of rotation of the planets, plus 180°. 


Determining the Radius Vector to the Sun 


After establishing the basic reference, and determining the 
vehicle latitude and longitude with respect to the sun, it is then 
necessary to obtain the distance of the vehicle from the sun. This 
can be accomplished using only one planet tracked with respect to 
the established reference system by a simple method proposed by 


Vertregct. en 


We can obtain from a memory drum or an ephemeris of 
heliocentric ecliptic positions of the selected planet as a function 
Of time; 

1. Longitude of the planet { Mes ) 


2. Latitude of the planet ( 5 ) ; 


5. Distance of planet from sun ( A ) 


To obtain the radial distance of the vehicle from the sun 


referring to Fig. 7,3, we can proceed as follows: 


ieee NCE 
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Los = ~A(E-O), (7-1) 


Fig. 7.2 Determination of heliocentric latitude and longitude 
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Fig. 7.3 Geometrical relationships determining radial 
distance from the sun 
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If the earth is one of the planets used for determining positon, 
and we use the plane of the ecliptic as the zero latitude reference 
plane, the solution is somewhat simplified due to the fact that 


earth heliocentric ecliptic latitude is at all times zero. (Fig.7.4 ) 
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Fig. 7.4 Geometrical relationships determining radial distance 
from the sun with earth as the reference planet 
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The determination of range of the vehicle from the sun can 
be made using only one planet after LA and AA are known. 
However, since we know appreciable errors exist in the values 
given for planetary positions as shown in Section 8.62, the use 


of two or more planets should tend to narrow the error down. 


Some Additional Methods for Position Determination 


A second method for angular determination of the sun and 
planets is photographic eH i‘). . .... By the taking of simultaneous 
pictures of a given small area (such as 5°) centered about the 
selected body, the position of the body with respect to the 
background of recognizable fixed stars can be obtained by taking 
measurements from the developed plate. This does not appear 
practical even for a manned vehicle, due to time lag, the need 
for continuous or very frequent data to effect smooth system 
performance and, in the case of a manned vehicle, an extra 


burden upon the crew. 


Similarly it has been suggested that the angular measurement 
of the diameter of solar bodies might be used for range determi- 
nation. This could be photoelectric or photographic in nature. 
Both these systems do not appear practical except possibly for . 
terminal guidance where the diameter of the body being approached 


becomes appreciable. 


Further postulated ranging methods are (1) measurement of 
incident solar radiation intensity, oe (2) measurement of reflected 


solar radiation from planets, and (3) radar ranging. 


Some basic arguments against use of the above methods for 


mid=-course ranging are as follows: 


1. The variation in intensity of solar radiation does not 
seem to be great enough for any but the most coarse range 
determination. Using Fig. 2.2 wesee that for a trip 
between Earth and Mars this variation in radiation is 
approximately i erg per square inch per second for a 


change in range of one mile. 


lid 


ie 


8 


2. Reflected solar radiation also has the same general 
gradiant as above in (1) plus the uncertain and possible 
widely varing amounts that might be reflected from a 

body do to changes in planetary atmospheric condition 


and orientation of lighted surface relative to the vehicle. 


5. Radar ranging has the same disadvantage as listed in 


Section 7.2 


Determination of the Shape of Actual Transfer Ellipse 


If we assume that after obtaining a series of acceptable 
vehicle positions we want to know the elements and orientation of 
the ellipse being traversed, we can proceed after first converting 
the differences in latitude and longitude of three selected positions 
into angular displacements along their great circle locus (@, - 6,). 


This can be accomplished by standard coordinate transformation. 


Determination of three of the basic elements of the unpowered 
orbit, 9, e, anda, can be accomplished by the use of three se- 
lected positions of the vehicle with respect to the sun. For any 
reasonable degree of accuracy the time interval between selected 
positions may have to be several days. A method of mathematical 
reduction which can be easily programmed for a computer has been 
suggested by Vertregt?! and is adapted herein to the nomenclature 


of this paper in Derivation Summary 7.1. 


Two assumptions are made in this method which are common 
to many proposed non-powered interplanetary mid-course tra- 


fecvories, | | hese are: 
1. Vehicle is following an unperturbed Keplerian ellipse, 


2. Sun is only gravitational force considered, and is one 


of the foci of this ellipse. 


From Derivation Summary 7.1 we obtain 


DETERMINATION OF THREE BASIC ELEMENTS OF ACTUAL 
ELLIPTICAL TRAJECTORY FROM OBSERVED DATA 


The following method of determining 9, e, and a is adapted 
from reference 67 to the nomenclature of this paper. 
Using a basic equation of any ellipse and referring to Fig.7.5 
we obtain the following: 
a (1 - e*) 
(5) 


R = 
A, l+ e cos 0, 


a (l- at 
Ane l1+ecos 0, (6) 
a(l- 2 
(7) 


R — 
A, l+ecos 04 


Combining (5) and (6) and (5) and (7) we get 


1 | 
Ra +R, e cos 6 = R, +R, e cos 05 (8) 
i I p 2 
| 1 
R, + Ry e cos a= R, +R, e cos 0. (9) 
l i 3 3 
Rearranging (8) and (9) 
1 1 
R, -R, = e(R, cos 6,-R, cos 6, ) 
A, A) A) i As 2 
“A, oo 
e = (10) 
- R cos O65 


DERIVATION SUMMARY 7.1 
(Page 1 of 3) 


117 


2oeds ul asdiija yayoor aut oct) SSI 


(L) 





NOMISHIY3d 
LIGdo 


JNOIHAA Wa 


Sw cre ee ee 





> 





N Hvd STOIHaA 7 NO! tiv 
JIDIH3A 4O ISO NOILOSrOud JIOIHBA 


NOILOW 4O NOILOSYIA 


H1LVd 
WLIGYO AFIDIHSA 


His 


— ? 
Ry -R, =e (R, cos 6,- R cos 6, ) 


3 1 1 43 : 
R - R R - R 
_ 43; 4 ; Bo 4, 
1 Wee 1 1 
Ry cos@)-R, cos ©, RA cos 6)- RK, cos 05 
] 3 I Ds 
(11) 


From (11) we obtain 


1 1 
RA (RA -R, ) cos 6,-R (R, - RA, ) cos @ 


1 As 2 4, 4, 3 ‘ 
1 

+R, (RA sek ) cos O. = 0 (12) 
3 4 2 
1 1 1 1 

0, = 0, + (8, - 6) (13) 
ee 

6, 7° 0, + (6,- 6,) (14) 


t ’ 1 1 

(0, - 6,) and (9, - 6,) have been determined along with R, , 
2 Hl 3 I 4, 

RA : R, by establishing the three instantaneous positions of the 


vehicle. ie can now substitute (13) and (14) into (12) 


RA a a cos 6, - Ra (RA - RA 


eos fo, +0, - a] 


3 
A, Sgy so Se fa, +65 - 6, 0 (15) 


Expanding the summed angle terms in (15) leads to 
1 


| 1 1 
R, (R, -R, )cos @0,-R, (R, -R, )Icos @, cos (6.-8@,)- 
4, As As ] A, ay As ] ie | 
t t t 1 | 1 
- sin a sin (6, - a) of Pe ae 37. cos @ cos (6, - 6) 
| 1 1 
- sin 6, sin (6. - a) = 0 (16) 


DERIVATION SUMMARY 7.1 
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t 
Solving (16) for 6, we get 


| q $ 8 
-R, )cos (@, -6,)+R, (R,-R, ) cos (6, - 8 
A, 2 i A, + A, 3 l 


R, (R -R )-R (R 
Sj A 4g Ag a, 4 
8, =x tan 


1 
MAG (RA - Ra ) sin (8, - @,)+ Ra aa an ) sin (6, “ 6,) 
3 as 2 a7) 





From (11), substituting the value obtained in (17) we obtain 


(18) 





Substituting the value obtained in (18) into (7) and rearranging 


leads to 


(19) 
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’ e | } 
R, (R -R, )-R, (R, -R, ) cos (@, -@,)+R, (R, -R, ) cos (6, - @ 
ao | A awe A An Naa Bo IN 3 
ik 


(R, - 
ay 


-R ; b) b : ’ i] 
(R >] Rh.) sin (6. = 6 )+R R, ) sin (0, - @,) 
bo A ae i A, ead 


(7-6) 


(7-7) 





Determination of Orientation of Actual Transfer Ellipse 


Two additional elements of the transfer ellipse, Q“ and i, are 
required for determining the orientation of the vehicle trajectory 
with respect to the reference plane. These can be obtained by 
applying the method suggested in reference oi using two positions 
of the vehicle. A rather obvious observation can be made that the 
further these two positions are from each other the more accurate 
will be the values of 22 and i obtained. However, since corrective 
action is more economical if performed early in the flight, it is 
apparent that £22 and i should be determined as early as possible. 

If the trajectory is satisfactory at this time more accurate values 


can be obtained at the next determination. 


From the first vehicle position, as illustrated in Fig.7.6 we 
can obtain the following: 


S, 1s to NN 


] 
O 
A, _ 2 05 =e) 
s/s 3} ) 1 
Z. ae sin Li 
] 
= = \ ae ( OO, _ 
S) = a cos AR, S) © Ra cos (90 n,) Ra, sin 1, 


12] 


asdij[a 1ajsuesz} JO udT}eWSITIOQ Oy “SIA 


(L ) 





ee 


therefore 


sin =Zy 
sin i = ————- (7-8) 
sin 1, 


For the second position we can likewise determine that 


sin[ n,t (0, - 6,) 


f 


Combining the above 


sin i : sin Sis, sin a, 


aaa sin[n, + (6,*- 6,)] sin n,cos (0, = O)+ cos 7, sin(6,- 0, ) 
Multiplying by sin ny and regrouping leads to 


sin LA Sin (@, = @,) 


1 ; z 


n) = tan — (7-9) 





& - / } j 
sin Ma - sin Ma cos (85 - 8.) 


n, can now be converted to the plane of the ecliptic. If we then 


subtract it from A A we have {. 


After determining n, we can then solve for the inclination of 


the orbit of the vehicle using equation (7-8) 


(7-10) 
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B, IN-FLIGHT CORRECTIONS 


7.10 Introduction 


It is considered inevitable that the rocket may fail to follow 
its prescribed path in space. Errors in initial conditions or 
unknown gravitational influences will cause deviations from the 


predicted trajectory. 


As soon as differences between actual position(as described 
in Part A above and predicted position become known, 
anew Situation arises. The decision must be made to either 
return to the previously calculated orbit by means of corrective 
thrust, or to recompute a new path which will be within the 
capabilities of the rocket. Time is an important consideration in 
this discussion Since corrective action must place the rocket on 
a path which will result in almost simultaneous coincidence of 


Mars and rocket. This is discussed further below. 


It 1s presumed that almost continuous position finding data 
will be available so that comparisons can, infact, be made 
between actual and predicted position. It is also reasonable to 
assume that deviations, hence corrections, will be incremental 
in nature. For the sake of simplicity, let us also assume that 
computation of a new trajectory, when required, can take place 
instantaneously. Under actual conditions, solution time can be 


allowed for with no great difficulty. 


The desirability of placing the rocket on a new path rather 
than forcing it to adhere to a previously defined path has been 
mentioned. The latter method implies. in effect, the application 


of continuous corrective thrust. Reference 26 has indicated an 


— 
Nh 
ey 


ml 


analytical approach to the thrust requirements necessary to fly 
a prescribed path. Until such time as fuel supply poses no 

problem, however, continuous corrective thrust is prohibitive. 
As little application of corrective thrust as possible is desired, 


since there is little available for expenditure. 


The objective of corrective action is notnecessarily to 
achieve a pinpoint interception of Mars, but rather to reach a posi- 
tion favorable for entry into orbit about Mars. This might 
conceivably be a position at some distance thousands of miles 
from the planet as well as on the planet itself, since the 
gravitational influence of the arrival planet and decelerating 
thrust will combine to make entry into orbit a less significant 
consideration. Therefore, some error can be tolerated. However, 
when and if deviations occur which are outside tolerance limit, 
corrective thrust becomes necessary. Obviously, the tolerable 


deviation is smaller at the early stages of flight. 


For purposes of analysis let us assume that the allowable 
tolerance has been exceeded and that corrective action is called 
for. The rocket is moving primarily under the attraction of a 
single body, the sun, which is justifiable in accordance with the 
calculations of section 3.6. We will also consider Mars as a 


massless point to be "intercepted" (in the sense defined above) 


at a specified time. 


The subsequent analysis parallels that of Lawden— 


although different in many respects. 


y 


Analytical weveropment 


er to , 8, define present position in space at time th . Let 
Ty. 0, Aetine ie position the See A would occupy at time t. had 


flown we correct path. Ty» 0, define the position of Mars = time 


126 


to, the planned time of interception. It is desired to intercept 
Mars at time to via a route which obviates the necessity of re- 
turning to the original trajectory. Figure 7.7 is helpful in the 


interpretation to follow. The route is a conic passing through 
evand: J, 
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/ f : es? swe 
Mars eg, E Ss get ce 
f @ co 
Trajectory Ones 
os oo 
Xs 
W 
ce ee els. | ae 4 reference line 
perihelion of predicted 


path 


Fig. 7.7 Interception at Originally Predicted 
Time, to 


Clearly, it is a restriction to be forced to intercept the 
planerat time, to, but results in a great simplification in the 


computer requirements. Let us pursue this analysis first. 


The conic path to be flown is presumed to be a segment of 
an: ellipse (because of velocity requirements) with sun as focus. 
The line of apsides of this new ellipse, whose elements are 
designated (aj, €5) is aligned at an angle w with the reference 


line of the originally prescribed ellipse (elements a1,€,). 


From the fundamental elliptical equation 4-1] 


ao = ey”) 
ie = ——\" = ak | nel (7-11) 
1+ e, cos ug - w) 


Ao (1 - a 
{> Se (7-12) 
1 +e, (6, - w) 
t t 
where bea To, 6 and 05 are known. 


LS 


Position in an elliptical orbit at any time t can be found from 


_ AT 
the relationship 


i e. + cos @ 
27t e 1-e- sin 0) ein I + 
: 7 


— 


rls 


I 
f 
toe cos 6, 1+ e cos 0, 


] 


where t is measured from t = 0 which occurs at perihelion passage 


of ellipse being considered... 


4. 
T, the period of the ellipse,can be expreeeede as a function 


of a, and is usually known. 
3/2 
(7-13) 


From these relationships, 


2 : ; 
Nil -e sin (0, - Ww) 








t 
2m to ey 9 1p &g + C0884 -w) T 
T 2 aS a =" lf) =". : ne ae 
2 1 +e, cos (6, -w) te, cos(6,-w) 
and ' t 
20 th egN 1 = an sin (0. = Wy) Fl eo + cos (Cae w) 7 
T ns Sil )| «6 aa Bay 
2 1 + ey cos (6. = ) 1+ e, cos (6, =i) 


where to ane ce are unknowns, measured from the reference line 


inclined at an unknown angle w to a known reference line. 


If the time cf travel from ae to l5 is constrained to be the 
originally planned interval of time from ry to Lo, the interval 
to ~ ty becomes a known quantity, and to and ty need not be known 


individually. 





! ' 
To @5N1 = en” sin (6. - w) €5 l-e,” (0 - w) 
SS ees ee 
O 7 
1+ €, cos (0, - w) 1+ e, cos (9 = w) 
? ? 
-l e, + cos (6, ~ q)) ie +cos(@.- Ww) 
3) SS SSS eee 
1+ € COS (65 -w) 1+ C5 cos(@.-w) 
(7-14) 


In summary, position finding data will provide lo OO: The 
! 
quantities ry» o, and lo. 
time. Thus, At, the time till interception, is known for any r 
f 


? 
G, . Equations 7 ll 7~!4involve four unknowns Ay, Co, W, and To, 


and are solvable by iterative methods. 


i] 
95 are known in advance as functions of 


piNce COrrections are assumed to be nerementalaenaraiee 
when compared with the elements of the path being flown(a,, ey 
Wo): a5, €y and w, may be replaced by a+ 6a, e + 6e, and 
w + dw. It is assumed Ty has been replaced by its value in 


terms of Ay: 


If these substitutions are made and small angle assumptions 
are valid (sin 6w = 6w, cos 6w= 1) equations in the three unknown 
6a, de, and dw result. ay? &5 and Wo (assumed zero, initially) 
can be found as described in the navigational section and are good 


first approximations toa, e andw. 


Successive improvements of the accuracy of 6a, 6e, and é6w 


can be made if desired. 


Corrective Acticn 


Since changes must inevitably be reflected in an alteration of 
velocity and angle of fire, we may again write two basic re- 


lationships 


(vio + ISEe Gigeee . i) 
O Lo a5 


b et 
Ww 
CO 


J 
e, sin (@, - w.) 
tan (¢. + 46.9 ——————— 
- © l1-e, cos (0, = wo) 
2 2 Z 
Replacing Ap» 5 and Wo by a+ éa, C5 +6e and w+ 6w, we 


arrive at 


ay 
2 
z ] 
=u -[om | tT, atéa I-¥< 
where wes is the velocity at the time corrective action is deemed 
necessary, and % is the corresponding angle of orientation as 


defined earlier. Similarly, direction must be altered to effect 


Ad, which can be written 


7 
_1 (e +6e) sin (6, - w ~ du) 
iy se 
O 


| 
1-(e+de) cos (65 -w dw ) . 


at+déa, e+d6e, and w + é6w are presumed to be known accurately 


as described in the preceding paragraphs. 


Corrective thrust must be actuated to produce the changes 


required. 


Interception at Newly Calculated Time 


Let us now consider the prospect of intercepting Mars at some 


time other than the originally planned time. 
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Fig. 7.8 Interception at Newly Calculated 
ime 


It is recognized at vA that present position, vee oe is un- 
Satisfactory,i.e. a continuation along present track will not 
achieve interception. Moreover, assume that rocket corrective 
thrust capabilities will not permit an interception at ro, Oy. the 
originally planned arrival point. We desire a conic (elliptical) 
path from lo: Q, to some arbitrarily selected Mars position, Ls, 
O'. which has been chosen becauSe it is within the rocket's thrust 
capabilities. This new position may also, in fact, represent the 
position of Mars which can be most economically reached. The 
choice of this position does present a formidable problem, since 
the optimization process is necessarily involved. The selection 
of a Mars position is not discussed here although its importance 


is emphasized. 


Presumably, there is on board a means of establishing Mars 
position as a function of some base time - a Mars Ephemeris of 
sorts. Fixing the desired interception point will fix the time, ta, 
of interception. ty is known from an accurate time piece. Thus, 


the interval of corrective travel becomes tt. - us . These times 


5 
are both based upon a common reference. The time interval is 


therefore determinable. 


Employing the previously developed reasoning procedures, 


2 

: a,(1 ~ 5 ) 

ae - : 
lt+e, cos (6 +w) 

2 

; ao(l - €5 ) 
a a 
l+e, cos (6.,+ Ww) 


2 ee 2) Sa 
T, ey l-e, sin (64+ 0) ep l-e, sin (6. + w ) 


At=t,-t => 

3 Or 27 : ; 
1+ C5 cos (6, +w) 1+ Co cos (0. + w) 
1 | 
lp) @g t cos (6, + w) 11 &gt cos (O6,+) 
- sin + sin oe 
1+ ey cos (0, +w) 1+ Co cos (0. +w) 
and 
20 a?! ; 
TT ae 


By the methods previously described, the elements of the 


new patha can be found and corresponding 6a, ée, and 


ey 2 
6w's determined. It is now a relatively simple matter to 
calculate the velocity and angular changes necessary to fly the 


corrected path. 


7,14 Summary 


It is seen that both approaches are essentially the same. 
From the standpoint of instrumentation, however, interception 
at time other than planned requires on board a knowledge of Mars 
position at all times and necessitates the selection of a new inter- 
ception point. The choice of this point is essentially an 


optimization process. 





Glaleu slide 6: 


GUIDANCE AND CONTROL 


8.1 General 


The instrumentation of any guidance and control system for 
interplanetary travel is extremely complex. As stated previ- 
ously, this paper studies the so-called "combination" trajectory 
(i.e. ~ an initial thrust period followed by unpowered free-fall 
flight; modified, if neaessary, by corrective thrusts). Primary 
emphasis has been placed upon what has been defined as the ''mid- 
course’ of this type trajectory. Although considerable thought 
has been devoted to the use of departure and arrival orbits, the 
applicability of this paper is not limited solely to interorbital 
trajectories. Whether the vehicle orbits prior to commencing 
mid-course travel or simply continues directly from the ascent 
from the surface of the Earth does not affect requirements of the 
vehicle velocity vector orientation and magnitude and other para- 


meters for the actual transfer period. 


This Chapter outlines in simple form a suggested Mid- 
course Guidance and Control system based upon the material 
developed in previous portions of this paper. The major problem 
of reducing the system to the least possible weight and size has 
been ignored except to the extent that simplification has been 
attempted where it will not detrimentally affect the overall ac- 
complishment of the mission (i.e. - the arrival of the rocket at 


Mars within certain limits). 


Time-sharing of a miniaturized low-power digital computer, 
a common ephemeris memory, and one basic time generating 


mechanism have been assumed. This appears justifiable since 


the elasped time between computations need not be of the order of 
fractions of a second, but can be several minutes or possibly hours. 
Digital computation is necessary due to the extremely low rates of 
change of rocket position, velocity, and acceleration; and the need 


for a high degree of accuracy. i 


The problems of instrumentation have been divided into 


several sections as follows: 


(1) The Basic System 
(2) The Tracking Section 
(3) The Computing Section 
(a) Tracking Line Computation 
(b) Prediction Computation 
(c) Comparing Function 
(d) Correction Computation 
(4) Thrust and Stabilization Section 
(5) Errors and Effects 


These sections are treated below. 


8.2 The Basic System 


The functional responsibilities of Guidance and Control have 
been segmented into three basic categories, as shown in Fig. 8-l. 


A brief general description follows: 


‘he Tracking Section - Maintains Lines-of-Sight upon 
prescribed bodies and provides 
Tracking Line data therefrom to 


the Computing Section. 


The Computing Section - Performs the necessary computations 
to the Tracking Lines to obtain 
Signals proportional to the cor- 
rections required to the rocket 
velocity vector to insure ac- 


complishment of the mission. 
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ihe Thrust and Stabi- 
lization Section - Provides vehicle orientation and 
corrective thrust in accordance 


with the received signals. 
® 


The aforementioned system sections will be examined in 
more detail in the following sections of this paper. The Computing 
section of the rocket involves such a complex problem that it has 
been subdivided into four parts, each of which will be discussed in 


its own Section of this Chapter. 


The Tracking Section 


The tracking section envisioned determines tracking lines 
from the vehicle to the Sun and planets measured in an helio~= 
centric ecliptic coordinate system. These angular measurements 
are then sent as electrical signals (TL) to the tracking line com-= 
puter in the tracking section. Two systems are presented which 
differ in method of space reference frame determination, number 
of star and planet trackers required, and in probable degrees of 


accuracy and reliability. 


System One 


The first system obtains space reference from star trackers 
rigidly mounted to a control member platform isolated from vehi- 
cular motion by gimbals. This platform is "fixed" in space and is 
maintained in this attitude by tracking error Signals generated in 
the fixed star tracking units, which are used to actuate the gimbal 
drive motors. A gyro package is mounted on the controlled member 
solely to assist in holding track during periods of radical maneuver 
and to aid in initial lock on. Three additional star trackers are 
mounted on the controlled member, each with its own tracking 
member drive system, as illustrated in Fig. 8.2 and 8.3. One 
tracking head determines the rocket-sun tracking line while the 


others determine two different rocket-planet tracking lines. 


famed 
cA 
ao 





SOLAR SYSTEM 
TRACKERS 






Se 
TRACKER VERTICAL PLANE ——% —— 


a Prey Ani ne 
7 is — - DRIVE a | 
>t ST 
‘y am - = 
7 Nom <4 y 


= | ce 
MIDDLE GIMBAL 


DRIVE MOTOR REFERENCE PLANE 


CONTROLLED MEMBER 


SPACE REFERENCE 
FIXED STAR TRACKERS 


OUTER GIMBAL S INNER GIMBAL GYROS 


DRIVE MOTOR DRIVE MOTOR 


BASE 
& 


Fig, 8.2 Illustrative schematic of tracking system 


with gyro assist 
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32 


, 33 


Angular measurements of these lines are made in heliocentric 
ecliptic coordinates, (see Fig. 7.1 ) and transmitted to two 
separate range computers within the tracking line computer of 


the computing section. 


system Two 


The second tracking system (Figs. 8.4 and 8.5) consists of 
a Space reference gyro package, and only three star trackers. 
The reduction in the number of trackers used is made solely to 
simplify construction at the expense of system accuracy Three 
orthogonally positioned single degree of freedom gyros are as- 
sisted by a single ‘fixed star’ tracker rigidly mounted to the 
controlled member. Smoothed error signals from this tracker 
are combined with the signals from the gyros. This resultant 
Signal is then used to actuate the gimbal drive motors maintaining 
the controlled member reference plane in proper orientation. 
Two Star trackers are then mounted on the controlled member, 
each with its own tracking member drive system. As in the first 
system. one tracker determines the rocket-sun line, while the 
second tracker determines a rocket-planet line. Angular measure- 
ments of these lines feed into a range computation section of the 


tracking line computer. 


Comparison of Systems One and Two 


It is obvious from comparing Figs. 8.2, 8.3, 8.4 and 8.9, 
that the second rnethod discussed in 8.32 is much Simpler to 
design and construct since only three stellar tracking heads are 
needed while the system discussed in 6. 3] requires six tracking 
heads. (Possibly this could be reduced to five if only two heads 
were used for space reference determination). The main draw 
back to system two is the dependence upon gyros for space refer- 
ence. These units are not, at present, Staple enouch No eapae. 
dictable to the degree required, for a time of flight of the length 


required for an JXarth-Mars transit. Section 6.64 discusses this 
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Fig. 8.4 Illustrative schematic diagram of tracking system 


with stellar assisted gyro space reference 
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8.4 
8.41 


particular aspect. 


A second drawback to the simplified system of section 8. 42 
is the use of only one planet The effect of this is felt in com- 


putation of O81 Sys! are and will be treated in,section 8. 41, 


The Computing Section 
Tracking Line Computation 

The tracking line computation section of the Guidance System 
receives the angular measurements of the tracking lines de-~- 
termined in the tracking section. Using these measurements and 
data from the Ephimeris Storage and Time Generator within the 
computing section, (RO) obs and Woy ane are determined. In 
addition the various separate elements of the ellipse are com- 
puted for complete identification of the actual trajectory being 


followed by the vehicle. 


The method of computation is geometric using discrete 
measurements, and not the integration of velocities and acceler- 


i, (2 The reasons for this 


ations experienced by the vehicle. 
are shown in section 8.6. Chapter 7 developed the equations 
necessary for the solution of this problem and illustrated the 


various geometric relations in Figs. 7. 2 through 7.6. 


The first tracking system discussed in section 8.3] has as 
outputs TL : The and TL: These are inputs to two separate 
range computers which determine two values of ROA using the 
equation of section’?.6., These two values are then averaged and 
smoothed for greater accuracy and then combined with the vehicular 
heliocentric latitude and longitude (LA: Ay) as presented by Tho . 
Since a series of RO, must be used in computation of the elements 
of the elliptical trajectory, these must be held ina memory storage 
unit, and fed to the trajectory elements computer upon demand. 


The determination of (Vo a) is made by a simple integration 


obs 


process and along with (R@),) is a direct output of the tracking 


obs 
line computer. The output of the trajectory elements computer 


are 9, a,e, andi observed. The simplified information flow 


poe 
aa 
ay 


42 


diagram of the above system is shown in Fig. 8.6. 


The second system is the same as that described above but 
uses only one range computer with the resultant loss of a degree 
of accuracy in the solution, although a system simplification is 
achieved which may be permissible. A simplified information 


flow diagram of this is illustrated in Fig. 8.7 


The Prediction Computation Section 


The Prediction Computation Section consists primarily ofa 
number of components collectively termed the Prediction Com- 
puter, (see Fig. 8.8). The purpose of the Prediction Computer 
is to continuously establish, as a function of time after te the 
position and velocity of the rocket relative to the reference body. 
The position and velocity so established may be termed "'pre- 
dicted" in the sense that they do not stem directly from obser - 
vations, but rather are values which should actually occur if no 
uncompensated events transpire to alter the rocket motion. The 
computer must establish these values without reference to any 
source of information from outside the rocket itself. This in no 
way prohibits the computer from sharing storage data, power 
Supply, time generation, and the like. with other components. 

If no corrections to the trajectory were envisaged all the outputs 


could be computed prior to launch, stored ina memory, and 


sampled as needed. Since corrections are envisaged the computer 


must possess the capability of solving values for the rocket po- 


sition and velocity during the entire time of flight. 


selecting the Sun as the center of our inertial reference 


frame we may write Equation 3-ll, 


Roa | 7 far t Goa + (Goa - Gag) t (Goa - Ge) + (Gor, - Se) 


(8-1) 
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om 


43 


mrom Big 2,8 “we determine 


-_ 


GQO < 107! Ph gees 


oS 


GQ@ < 107? Pleas (8 -2) 


GIO = 1078 BE Gee" 


Hence the equation of motion to be solved continuously by the 


Prediction Computer is 


Roa | = tip + Goa, t Coa + Gqa t+ Gg, (8-3) 


The threshold sensitivities of the components will determine 
the Phases referred to in Section 3.6. Below these thresholds 
a particular component will be inactive and the problem will be 


reduced to the appropriate one of Equations 3-14, 3-15, or 3-16. 


Referring to Fig 38 6 the functions of the sub-~components 
are almost seli-explanatory when examined with regard to 
Equation 8-3. The Summer forms the quantity representing the 
acceleration of the rocket with respect iomiewoun mene acceler. 
ation is integrated twice to form the rocket position vector, 
which is fed-back to obtain values for those gravity fields con- 
sidered in the equation. Thrust and rocket mass form the 


remaining inputs tothe Summer . 


Comparator Section 
ee Emenee 


The comparator section may be considered an integral part 
of the main computer’ Its function is to compare observed 

: : (Q . + . 
position in space, RO) obs: emanating from the tracking 
cormputer, wiih predicted position in space, (RO, PRE derived 


from the prediction computer. Deviations of observed position 


14g 


magnitude, form the basis for corrective action. The test of 
"sufficient magnitude’ is a function of the position of the vehicle 
with respect to Mars. Obviously, the permissible deviations 
would be small during the early stages of flight, when the position 
of the rocket relative to Mars is quite larpe. Similarly, as this 
relative position becomes less, tolerances become larger and 
errors assume a lesser importance. A deviation computer 
weighs or assesses the rocket position relative to Mars and 
creates this standard of measurement. Fig. 8.9 is a functional 


representation of this section. 


If actual deviations in position, (D,) RO,: exceed those 
permissible deviations, (D.) Rp; the correction comparator 
is assumed to emit a signal actuating the corrector section. 
This will automatically cause initiation of the computation of 


the new trajectory. 


8.44 The Corrector Section 


The corrector section can properly be considered a most 
complex and critical part of the computer. Upon its proper 
performance will depend an accurate approach to the objective 


planet or the alternative, a body in free orbit in space forever. 


The corrector section is continually supplied with observed 
rocket position, (ROA) obs: and the elements of the observed 
trajectory, as well as the position of Mars with respect to the 
Sun as afunction of time, R@At). These quantities form the 
basis for the computation of a new path to the sought-after planet. 


A discussion of this general problem has appeared in Chapter 7. 


The initiation of operations within this section will be de- 
pendent upon receipt of the comparator section actuating signal, 
Se(Act). 


The modified trajectory is asSumed chosen through an 


optimization process which minimizes fuel expenditure, time, 
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or some predesigned parameter. Once this selection process 

is completed, signals proportional to the required velocity 
changes can be converted into control system commands, ST, 
which will eventually translate the changes into corrective 


action. (See Fig. 8.10.) 


From a practical viewpoint, however, the selection of an 
alternate path in flight can be expected to pose a difficult 
problem. Size and weight considerations will necessitate an 
absolute economy in functional employment of components. ven 
then, the optimization problem is complex and may well repre- 


sent a major obstacle. 


.5 Thrust and Stabilization 


The thrust and stabilization unit is essentially a modification 
of the familiar control loop, altered because of operating en- 
vironment to include rocket thrust device rather than control 
surfaces. A signal proportional to the required velocity actuates 
the thrust device which is positioned prior to firing. Several 
means of directional contro! are available for corrective thrust 


requirements: 


(1) Movable vanes in the jet stream 
(2) Small, individual variable thrust jets 
(3) Universally mounted main jet engine 


(4) Gimbal mounted rocket thrust chamber 


Positioning of the thrust device would be in response to the 
correct or section output signal, as would the controlled 
metering of propellant to insure a velocity change of proper 


magnitude. 


The accelerated motion of the rocket produced by the thrust 
device would then be sensed by an accelerometer which com- 


pares. in effect, commanded thrust, with resulting thrust. 


The stabilization of the rocket in space is another problem 


lod 


of no little difficulty. Some of the proposals for rocket orien- 
tation have included micro-jet control, gyroscopic torque 
stabilization, and flywheels. Many authors combine vehicle 
motion and the stabilization although they distinctly differ and 


° . c &® 
require functional separation. 


Fig. 8.11 is a line functional information flow diagram of a 


proposed thrust and stabilization sections. 


8.6 Errors and Effects 


It has been stated previously that the unpowered elliptical 
orbit trajectory for transfer between planets is subject to large 
errors in terminal position caused by very small perturbing 
forces during transit and also very slight errors in velocity 
vector orientation at the end of the initial thrust period. If, as 
has been assumed in this paper, corrective action must be taken 
during mid~course, the reduction in terminal error will be a 
function of time remaining till intercept, accuracy in determining 
the actual trajectory being traveled and the accuracy of determi- 
nation and orientation of the corrective velocity vector. We will 
now show the magnitude of some of the terminal errors caused 


by various assumptions and inaccuracies. 


8.61The Assumption of Circular Coplanar Orbits 


The assumption of coplanar circular orbits has been made 
by numerous authors for the sake of simplicity in presenting 


basic ideag 22° 202 31 40; 62, 68 


and others. If we compare the 
orbits of Earth and Mars and transfer between these orbits we 


find the following to be true for the worst possible cases. 


At "favorable opposition" (once every 16 years) Mars is 
Peeaieeieetron Earth. At ‘unfavorable opposition" Mars is 
63 x ik miles from Earth. This gives a maximum actual vari- 
ation in distance between orbits of 28 x 10° miles (.3 AU), while 


the so called average distance is only 49 x 10° miles (.523 AU). 
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This is a 57% variation which is ignored with the assumption of 


circular orbits. Using the coplanar assumption (not taking the 





actual orbital inclination of 1° 51'' into account) we introduce a 
terminal vectical error of 4.57 x 10° miles for the worst possible 
&® 


Case. 


8.62Errors in Astronomical Distance 


The basic unit of astronomical distance is the Astonomical 
Unit (AU). This is determined by the mean range of the sun from 
the Earth and by international agreement (Paris, 1911) has been 
established to be 93 x 10° miles. However, this distance is 
based on angular measurements of the diameter of the sun. Due 
to limits of resolution, limited base lengths for triangulation 
and other effects, this basic unit is thought to have a probable 


error of 10,000 miles. 5 


This uncertainty is therefore 
carried over to the measurement of positions of solar bodies 

and hence to vehicular orientation in solar space. Until more 
accurate astronomical measurements are available a rather large 
residual error in spatial position will have to be accepted and 


worked with. 


8.63Errors in Initial Velocity Vector 


If we now take a typical trajectory from Earth to Mars, 
we can see the effects of errors in initial velocity vector on 
terminal positions. Assume a trajectory of arc length 354 x 10° 
miles such as used in section 4.14 . This trajectory has a tirae 
of flight of 172 days, and an orbital inclination of 2.51°. An initial 
velocity error of 1 foot per second for this trajectory will give 
an error in arc length alone of 2820 miles. This is neglecting 
the effects on the actual elliptical properties of the orbit which 
increase the actual miss distance even more. Using an average 
orbital velocity of 25.8 miles per second, this would be equiva- 


lent to an error in arrival of 118.5 seconds, at which time Mars 


196 


(traveling approximately 15 miles per second) would be 1780 miles 


behind the initially proposed intercept point. 


An error in inclination or angular orientation of the initial 
Vel@ely vector of .0] degrees gives a terminal error nowmal 10 
the trajectory of 18,600 miles. We can also approach initial 
angular error from the standpoint of what will be the requirement 
for insuring a terminal error of not more than 1,000 miles. This 


will necessitate an accuracy within 1.94 second of arc. 


J.G. Porter*” discusses the effects of small initial velocity 
vector errors on the shape of a simple Earth-Mars elliptical 
trajectory. This trajectory starts at its perihelion and arrives 
at Mars at aphelion and has a basic time of flight of 281 days. Porter 
concludes that in pinata! Jee this trajectory to arrive at Mars 
within a distance of 50,000 miles, the velocity must be controlled 
to .001 miles per second or one part in 18,000. This would mean 
that to get within 1,000 miles of the desired terminal position the 


initial error would have to be less than .106 feet per second. 


For this same trajectory Porter concludes that an error 
of 12 seconds of arc in initial velocity vector orientation will lead 
to a 50,000 mile error at termination. This means that for a 
1,000 mile error we would have to be within .4 of a second of 


aera purn Out. 


The major conclusion which can be drawn from a look at 
errors in initial velocity vector orientation is that a simple 
unpowered trajectory with no corrective thrust applied during 
mid-course will not suffice for interplanetary travel. This was 
previously stated in section,4.4. However, of the two errors con- 
sidered, velocity and angular, the errors in angular orientation 


are the more critical. 


8.64 Mid--Course Errors 


We can now look at errors in some instruments which might 


5 ara 
i 7 
a & 


be used within the vehicle for guidance purposes such as star 
trackers, accelerometers, gyro packages, and associated 


equipment. 


If we presume the need for a maximum star tracker error 
equal to or less than 1, 000 miles in determining the positions of 
solar bodies 1 AU distant, we obtain a limit of 2.22 seconds or 
.0382 minutes of arc. This is a maximum allowable closed loop 
error composed of errors in reference plane, tracking head, 


and all associated equipment. 





Another approach to this problem of star tracker error is 
to assume that each sighting is a random event with Gaussian 
distribution of the error in angular measurement. If we now 
postulate that we have a "maximum" error of -5 seconds of arc 
in the values obtained from a star tracker and that 99.7% of the 
measurements will lie within these values, 95% will be in error 
less than 3.33 seconds while 68% will be in error less than 1.67 
seconds. Since an error of 1 second of arc is equivalent to 450 
miles at a distance of 1 AU, we can assume that 68% of the 
measurements would have an error less than 750 miles, and 95% 
would have an error less than 1, 500 miles when sighting bodies 
1 AU distant. It appears, then, if a number of readings are 
taken over a short period of time with proper smoothing and 
averaging, suitably accurate values can be obtained from tracker 
systems which may by themselves have a considerable possible 


Crror.. 


To determine an upper limit of the threshold of sensitivity 
to acceleration required for the trajectory of Section 4.14 of this 
paper we can again use a maximum error of 1,000 miles at 
termination. This corresponds to 4.75 x io fecinper sccond 
for a time of flight of 172 days. For this same trajectory the 
maximum allowable drift of a purely "inertial" + guidance system 
to insure vertical displacement error of not more than 1, 000 miles 


would be 1.64 x 10°? mr per day. This is determined using the 


trajectory arc length of 354 x 10° miles and time of flight of 172 
days. 


The primary conclusion that can be drawn from the above 
simplified mid-course guidance system errors is that although 
completely self contained inertial guidance systems are highly 
desirable, the required accuracies are rather high; particularly 
if return trips are contemplated without rezeroing the system. 
However, with the addition of star tracking units, suitably 
accurate guidance can be obtained which is not limited by or 


dependent upon time of flight. 
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CHAPTER 9 


CONCLUSIONS AND RECOMMENDATIONS 
FOR PURE ox 


Intreaguction 


The contents of this paper lead to several conclusions and 
recommendations of interest to those pursuing the problem of 
interplanetary travel. Many areas included in this thesis are 
presented with a view toward future use by students of the field 
and lead to a better definition of the problem. The bibliography 
represents a cornprehensive literature search of available un- 


classified related material. 


Conclusions 


(1) During the time of flight of a rocket from Earth to 
Mars the center of ithe Sun may be considered an 


inertial point. 


(2) During a major portion of the trajectory the problem 
may be considered one of two bodies, with some 


perturbations. 


(3) Threshold sensitivities and drift uncertainties of 
current gyroscopic mechanisms preclude the use of 
; ree ll ee 
a pure inertial navigation system at this time, par- 


ticularly when a return trip 1s contemplated. 


(4) During the time of flight of the rocket the motion of the 
"fixed'' stars is so small that they provide a means of 


establishing a space reference frame. 


(5) Determination of rocket present position in space within 


the limits proposed in this paper lies within the capa- 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


bilities of the current ''state of the art". 


Within the limits of the assumptions made, the 
"combination" trajectoryfie ,an initial thrust period 
followed by free fall flight; modified, if necessary, by 
corrective thrust)appears to offer the most immedi- 


ately feasible method of travel to Mars. 


Rocket velocity vector orientation appears to be more 


critical than magnitude. 


The Hohmann aliliseees provides the minimum required 
velocities to get to Mars but if rocket capabilities 
include V.__ (as defined within this paper) of 10, 000ft/sec 


or more a wide range of trajectories is available. 


The Rocket Plane Concept provides a means of reducing 
the trajectory to a two-dimensional problem, thereby 
greatly simplifying rocket construction requirements. 

It is emphasized that certain assumptions often made 

in the literature to reduce the problem to two-dimensions 
are unrealistic when viewed in the light of actual 
condition. The Rocket plane concept does not assume 
two-dimensions but selects the reference plane such 

that deviations from this reference plane are negligibly 
small. 


On a kinetic energy basis a certain optimum departure 
orbit exists for any given trajectory to Mars. Similarly 
an optimum arrival orbit exists about Mars. Additionally 
these optimum orbit radii are characterized by the 
property that for a relatively small increase in available 
velocity a satisfactory departure orbit radius may be 
used which is much nearer to the Earth surface than 


the optimum orbit radius. 


Instrumentation for travel between planets appears to 


lie within the capabilities of the present "state of theart"’. 





9.3 Recommendations for Further Study 


(1) Digital computer trajectory analysis should be 


instituted covering: 


(a) Required parameters for all trajectories whose 


initial rocket velocity lies within reasonable limits. 


(b) The effect of variations in initial conditions on 


linear miss distance at Mars. 


(c) Perturbative effects of solar system masses which 
were ignored in this paper due to their small 


magnitude. 


(2) Additional study should be instituted on the design and 
construction of free-fall space reference stabilization 


systems suitable for interplanetary travel. 


(3) Studies should be commenced which would lead to the 
establishment of optimization standards for trajectories 
to Mars on the basis of both fuel expenditure and 


harmeiic enercy ; 


(4) <A study of permissible deviations from prescribed path 
during flight allowing terminal accuracy within pre=- 


scribed limits should be instituted. 
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SYMBOLS 
DEFINITION SUMMARY 1 


symbol 


AiO DiC OS LP QOAO 


Definition 


sun 
Moon 
Earth 
Mars 
Jupiter 
saturn 
Uranus 
Neptune 
Pluto 
Mercury 
Venus 
star 


semi major axis of any ellipse 
semi minor axis of any ellipse 


The magnitude of the effective 
exhaust velocity vector measured 
with respect to the walls of the 


rocket nozzle 


Eccentricity of any ellipse: 


defined by the relation 
b2 
e =] a 
a 
Unit field forces - forces per 


a> 5 


unit mass whose existence is due 


to the presence of a field 


Symbol 


Definition 


Unit non field force - forces per 
unit mass whose existence is due 


to all causes other than fields 


m® 


Inclination of orbit - the angle 
measured from reference plane 

to plane of orbit, in a plane 
perpendicular to the line of nodes; 
generally positive where measured 


tH 


in a direction "above" the refer- 


ence plane 





A variable indicator which may 
take on identities as propounded 


within the text 
Instantaneous mass of rocket 
An arbitrary planet 


Magnitude of radius vector from 
the reference focus of an ellipse 
to the indicated point on that 
ellipse. Subscript j refers to 


indicated quantity 


Radius of the circular orbit about 
the indicated body 


~ 


The independent variable time 


Burning time - that period of 
time wherein a given amount of 


propellant is consumed 


The time at which an imaginary 
rocket would have arrived at Mars 
future position if it had departed 


aphelion of its ellipse at t = t 
O 
NGte: t, = tae eat 
16 e d 


op) 


symbol 


AU 


0 | 


| 


_—— 


~] 


Definition 


The time at which an imaginary 
rocket would have arrived at 
Barth present position if it had 
departed aphelion of its ellipse 
at t = ve 


Time of flight 
Time at "j" instant; subscript 


"j' refers to indicated quantity 


The instant of time at whicha 


given operation commences. 
Angle 


Astronomical Unit - the distance 
associated with the mean parallax 
of the sun. By international a- 
greement (Paris 1911), this was 

set as 8.80" of arc and corre- 
sponds to about 93,000,000 miles. 
In this text the numerical value of 
the AU will be taken as 93, 000, 000 


miles. 
Drag ~- as conventionally defined 


The ecliptic plane - the plane 
defined by the loci of points 
representing the instanianeous 
position of the earth throughout 
the year measured with respect 


to the sun 
Eccentric Anomaly 


Force - as conventionally defined 


Symbol 


Q| 


GM. 
J 


Imp 


KE 


Ses 1S 


Definition 


Gaussian gravitational constant 
approximately equal to: 


6.67 x10. em” /eram (m) sec” 


Directed acceleration due to 
gravity field of the indicated 
body. A second subscript refers 
to that body experiencing this 
directed acceleration, Subscript 


''j'' refers to indicated body 


The product of the Gaussian 

gravitational constant and the 

mass of the indicated body. Sub- 
tt.tt 


script j refers to indicated 
body 


Impulse - the integral of the 
magnitude of thrust over the 


time of its duration 


Kinetic Energy - as conventionally 
defined 


Celestial heliocentric latitude - 
angle measured in a plane perpen- 
dicular to the ecliptic plane from 
the ecliptic plane to the line drawn 
from the sun to the indicated point; 
measured positively toward ce- 


lestial north 
Lift - as conventionally defined 
Total mass of propellant 


Gross mass of rocket at instant 


Of Tie ie 
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Symbol Definition 


NN! Line of nodes - line drawn through 


the ascending and descending nodes 


PE Potential Energy - as conventionally 
defined 

PR Payload Ratio - PR=1- ¢ 

R Vector from the reference point 


to the indicated point 


R Magnitude of R 
fae Thrust - reaction force on the 
‘ rocket or vehicle unit structure 
At Total Energy - as conventionally 
defined 
Yo; Velocity vector of the indicated 


body with respect to a heliocentric 


non rotating coordinate system. 


Velocity vector of the rocket 
relative to the indicated body, 
and measured in a heliocentric 


non rotating coordinate system 


Velocity vector of the rocket in 
circular orbit about the indicated 
body measured in a body centered 


non rotating coordinate system 


Velocity vector of "escape'' from 
the gravity field of the indicated 
body, equal to 


N2 x V. 


js 
measured in a body centered non 


rotating coordinate system 
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Symbol 


VG A)g 


inc 


Definition 


Velocity vector of ine rocrer with 
respect to the indicated body 
taking into account the kinetic 
energy required to overcome the 
€irect of the gravity tielaor tne 
indicated body, measured in 
heliocentric non rotating coordi- 


nate system 


That increment of Vuale which 
must be provided by the rocket 


Angle measured from the x-axis 


to the thrust vector 


Angle measured from the y-axis 


to the thrust vector 


Angle measured from the z-axis 


to the thrust vector 


x,y, and z axes represent an 
orthogonal coordinate set. The 
thrust vector has its origin at 


the origin of this coordinate set 
M 
Propellant loading factor = aa 
O 
@-w for the orbit of indicated 
body 


The angle, measured in the rocket 
plane of motion, between a body- 
sun line and rocket ellipse aphelion- 
sun line; positive angles when di- 
rection of measurement is opposite 


to rotation of planets about the sun. 


symbol 


8 (contd) 


Definition 


Prime (') indicates measurement 
from perihelion vice aphelion with 
direction of measurement positive 
when same as rotation of planets 


about sun 
9 at any time j 


Gatt=ty 


G atrt St 
c 


6 of Mars att=t, 


Celestial heliocentric longitude 

- angle measured in the ecliptic 
plane from the reference direction 
Aries (7°), in the direction of 
rotation of planets about the sun, 
to the projection on the ecliptic 
plane of the line drawn from the 


sun to the indicated point. 


Longitude of perihelion of in- 
dicated orbit 


The angle, measured in the 
rocket plane of motion included 
between the velocity vector of 

the rocket or vehicle and a line 
perpendicular to the radius vector 


from the sun to the vehicle 


gat t=t, 


@ at t=t, 


The angle measured in the rocket 


plane of motion included between 


Symbol Definition 


Ye A (contd) V@q and the projection of Voge 
on the rocket plane. 


SA The angle measured in the rocket 
plane of motion, included between 
VA and the projection of Vod 


on the rocket plane. 


wd. The angle included between 
perihelion and the line of nodes 
measured in direction of ro- 


tation of planets for indicated 


orbit. 
A Rocket or vehicle 
ae The period - elapsed time between 


repetitive events 


WP The vernal Equinox (the first point 
of Aries) defined as the direction 
indicated when projection of the 
line of intersection of the earth 
equatorial plane and the earth 
ecliptic plane passes through 
the center of the sun; this action 
occuring as the sun apparently 
passes from south to north of the 


equator 


Aes Longitude of ascending node of 


indicated plane 


G; Longitude of descending node of 


indicated plane 


Ascending node - celestial heli- 


ocentric longitude occupied by a 


symbol 


* 
Loe | 
* ©« 
Need 


Definition 


body passing from ''below”’ to 
"above'' the ecliptic plane. 
Descending node - Ascending 


node plus 180° " 


Above’ ecliptic 
defined as that part of space 
wherein an observer would note 
counter -clockwise rotation of 


planet. 
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Fig. 4.12 Earth angular position measured from Earth aphelion : 
as a function of time (Overlay #1) 
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